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. Abstract 

Stochastic semiclassical gravity of the 90 's is a theory naturally evolved from semiclassi- 
cal gravity of the 80's and quantum field theory in curved spacetimes of the 70's. Whereas 
I semiclassical gravity is based on the semiclassical Einstein equation with sources given by the 

. expectation value of the stress-energy tensor of quantum fields, stochastic semiclassical gravity 

^ I is based on the Einstein-Langevin equation, which has in addition sources due to the noise 

kernel. The noise kernel is the vacuum expectation value of the (operator-valued) stress-energy 
bi-tensor which describes the fluctuations of quantum matter fields in curved spacetimes. In the 
bJQ. first part, we describe the fundamentals of this new theory via two approaches: the axiomatic 

and the functional. The axiomatic approach is useful to see the structure of the theory from the 
framework of semiclassical gravity, showing the link from the mean value of the stress-energy 
' tensor to their correlation functions. The functional approach uses the Feynman- Vernon influ- 

^ . ence functional and the Schwinger-Keldysh closed-time-path effective action methods which are 

convenient for computations. It also brings out the open systems concepts and the statistical 
and stochastic contents of the theory such as dissipation, fiuctuations, noise and decoherence. 
We then focus on the properties of the stress energy bi-tensor. We obtain a general expression 
for the noise kernel of a quantum field defined at two distinct points in an arbitrary curved 
spacetime as products of covariant derivatives of the quantum field's Green function, and show 
from this that the trace anomaly of the noise kernel is zero for massless conformal scalar fields. 
In the second part, we describe three applications of stochastic gravity theory. First, we con- 
sider metric perturbations in a Minkowski spacetime. We offer an analytical solution of the 
Einstein-Langevin equation and compute the two-point correlation functions for the linearized 
Einstein tensor and for the metric perturbations. Second, we discuss structure formation from 
the stochastic gravity viewpoint, which can go beyond the standard treatment by incorporat- 
ing the full quantum effect of the inffaton ffuctuations. Third, we discuss the backreaction of 
Hawking radiation in the gravitational background of a quasi-static black hole (enclosed in a 
box). We derive a ffuctuation-dissipation relation between the fiuctuations in the radiation and 
the dissipative dynamics of metric fluctuations. 
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1 Overview 



Stochastic semiclassical gravity is a theory developed in the Nineties using semiclassical gravity 
(quantum fields in classical spacetimes, solved self-consistently) as the starting point and aiming at 
a theory of quantum gravity as the goal. While semiclassical gravity is based on the semiclassical 
Einstein equation with the source given by the expectation value of the stress-energy tensor of 
quantum fields, stochastic gravity includes also its fiuctuations in a new stochastic semiclassical or 
the Einstein-Langevin equation. If the centerpiece in semiclassical gravity theory is the vacuum 
expectation value of the stress-energy tensor of a quantum field, and the central issues being 
how well the vacuum is defined and how the divergences can be controlled by regularization and 
renormalization, the centerpiece in stochastic semiclassical gravity theory is the stress-energy bi- 
tensor and its expectation value known as the noise kernel. The mathematical properties of this 
quantity and its physical content in relation to the behavior of fluctuations of quantum fields in 
curved spacetimes are the central issues of this new theory. How they induce metric fluctuations 
and seed the structures of the universe, how they affect the black hole horizons and the backreaction 
of Hawking radiance in black hole dynamics, including implications on trans-Planckian physics, are 
new horizons to explore. On the theoretical issues, stochastic gravity is the necessary foundation to 
investigate the validity of semiclassical gravity and the viability of inflationary cosmology based on 
the appearance and sustenance of a vacuum energy-dominated phase. It is also a useful beachhead 
supported by well-established low energy (sub-Planckian) physics to explore the connection with 
high energy (Planckian) physics in the realm of quantum gravity. 

In this review we summarize major work and results on this theory since 1998. It is in the nature 
of a progress report rather than a review. In fact we will have room only to discuss a handful of 
topics of basic importance. A review of ideas leading to stochastic gravity and further developments 
originating from it can be found in Refs. [1, 2]; a set of lectures which include a discussion of the 
issue of the validity of semiclassical gravity in Ref. [3]; a pedagogical introduction of stochastic 
gravity theory with a more complete treatment of backreaction problems in cosmology and black 
holes in Ref. [4]. A comprehensive formal description of the fundamentals is given in Refs. [5, 6] 
while that of the noise kernel in arbitrary spacetimes in Refs. [6, 7, 8]. Here we will try to mention 
all related work so the reader can at least trace out the parallel and sequential developments. The 
references at the end of each topic below are representative work where one can seek out further 
treatments. 

Stochastic gravity theory is built on three pillars: general relativity , quantum flelds and 
nonequilibrium statistical mechanics. The first two uphold semiclassical gravity, the last two span 
statistical field theory. Strictly speaking one can understand a great deal without appealing to 
statistical mechanics, and we will try to do so here. But concepts such as quantum open sys- 
tems [9, 10, 11] and techniques such as the infiuence functional [12, 13] (which is related to the 
closed-time-path effective action [14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]) were a great help in 
our understanding of the physical meaning of issues involved toward the construction of this new 
theory, foremost because quantum fluctuations and correlation have become the focus. Quantum 
statistical field theory and the statistical mechanics of quantum field theory [26, 27, 28, 29] also 
aided us in searching for the connection with quantum gravity through the retrieval of correlations 
and coherence. We show the scope of stochastic gravity as follows: 

^We will often use the shortened term stochastic gravity as there is no confusion as to the nature and source of 
stochasticity in gravity being induced from the quantum fields and not a priori from the classical spacetime. 
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I. Ingredients 

A. Prom General Relativity [30, 31] to Semiclassical Gravity. 

B. Quantum Field Theory in Curved Spacetimes [32, 33, 34, 35]: 

1. Stress-energy tensor: Regularization and renormalization. 

2. Self-consistent solution: Backrcaction problems [36, 37, 38, 39, 40, 41, 42]. 

3. Effective action: Closed time path, initial value formulation [14, 15, 16, 17, 18, 19, 20, 21, 

22, 23, 24, 25]. 

4. Equation of motion: Real and causal. 

C. Nonequilibrium Statistical Mechanics: 

1. Open quantum systems [9, 10, 11]. 

2. Influence Functional: Stochastic equations [12]. 

3. Noise and Decoherence: Quantum to classical transition [43, 44, 45, 46, 47, 48, 49, 50, 51, 

52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76]. 

D. Decoherence in Quantum Cosmology and Emergence of Classical Spacetimes [77, 78, 79, 80, 

81, 82, 83]. 

II. Theory 

A. Dissipation from Particle Creation [21, 22, 23, 24, 25, 84]. Backreaction as Fluctuation- 

Dissipation Relation (FDR) [85]. 

B. Noise from Fluctuations of Quantum Fields [1, 86, 87]. 

C. Einstein-Langevin Equations [87, 88, 85, 89, 90, 91, 92, 5, 6, 93]. 

D. Metric Fluctuations in Minkowski spacetime [94]. 

III. Issues 

A. Validity of Semiclassical Gravity [95, 96]. 

B. Viability of Vacuum Dominance and Inflationary Cosmology. 

C. Stress-Energy Bitensor and Noise Kernel: 

1. Regularization Reassessed [7, 8]. 

IV. Applications: Early Universe and Black Holes. 

A. Wave Propagation in Stochastic Geometry [97]. 

B. Black Hole Horizon Fluctuations: Spontaneous/ Active versus Induced/Passive [98, 99, 100, 101, 

102, 103, 104, 105, 8]. 
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C. Noise induced inflation [106]. 

D. Structure Formation [108, 109, 110, 111, 112]. Trace anomaly-driven inflation [113, 114, 115]. 

E. Black Hole Backreaction as FDR [116, 117, 118, 119, 120, 121, 122, 123]. 

V. Related Topics 

A. Metric Fluctuations and Trans-Planckian Problem [102, 103, 104, 105, 124]. 

B. Spacetime Foam [125, 126, 127, 128, 129]. 

C. Universal 'Metric Conductance' Fluctuations [130]. 

VI. Ideas 

A. General Relativity as Geometro-Hydrodynamics [131]. 

B. Semiclassical Gravity as Mesoscopic Physics [132]. 

C. From Stochastic to Quantum Gravity: 

1. Via Correlation hierarchy of interacting quantum fields [2, 27, 29, 133]. 

2. Possible relation to string theory and matrix theory. 

We list only the latest work in the respective topics above describing ongoing research. The 
reader should consult the references therein for earlier work and the background material. We do 
not seek a complete coverage here, but will discuss only the selected topics in theory, issues and 
applications. We use the (+, +, +) sign conventions of Refs. [30, 31], and units in which c = h = 1. 

2 From Semiclassical to Stochastic Gravity 

There are three main steps that lead to the recent development of stochastic gravity. The first step 
begins with quantum field theory in curved spacetime [134, 32, 33, 34, 35], which describes the 
behavior of quantum matter fields propagating in a specified (not dynamically determined by the 
quantum matter field as source) background gravitational field. In this theory the gravitational field 
is given by the classical spacetime metric determined from classical sources by the classical Einstein 
equations, and the quantum fields propagate as test fields in such a spacetime. An important process 
described by quantum field theory in curved spacetime is indeed particle creation from the vacuum, 
and effects of vacuum fluctuations and polarizations, in the early universe [135, 136, 137, 138, 139, 
140, 141, 142, 134, 143, 144], and Hawking radiation in black holes [145, 146, 147, 148, 149]. 

The second step in the description of the interaction of gravity with quantum fields is back- 
reaction, i.e., the effect of the quantum fields on the spacetime geometry. The source here is the 
expectation value of the stress-energy operator for the matter fields in some quantum state in the 
spacetime, a classical observable. However, since this object is quadratic in the field operators, 
which are only well defined as distributions on the spacetime, it involves ill defined quantities. It 
contains ultraviolet divergences the removal of which requires a renormalization procedure [134, 150, 
151]. The final expectation value of the stress-energy operator using a reasonable regularization 
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technique is essentially unique, modulo some terms which depend on the spacetime curvature and 
which are independent of the quantum state. This uniqueness was proved by Wald [152, 153] who 
investigated the criteria that a physically meaningful expectation value of the stress-energy tensor 
ought to satisfy. 

The theory obtained from a self-consistent solution of the geometry of the spacetime and the 
quantum field is known as semiclassical gravity. Incorporating the backreaction of the quantum 
matter field on the spacetime is thus the central task in semiclassical gravity. One assumes a general 
class of spacetime where the quantum fields live in and act on, and seek a solution which satisfies 
simultaneously the Einstein equation for the spacetime and the field equations for the quantum 
fields. The Einstein equation which has the expectation value of the stress-energy operator of the 
quantum matter field as the source is known as the semiclassical Einstein equation. Semiclassical 
gravity was first investigated in cosmological backreaction problems [36, 37, 38, 39, 40, 41, 42, 154, 
155, 156], an example is the damping of anisotropy in Bianchi universes by the backreaction of 
vacuum particle creation. Using the effect of quantum field processes such as particle creation to 
explain why the universe is so isotropic at the present was investigated in the context of chaotic 
cosmology [157, 158, 159] in the late seventies prior to the inflationary cosmology proposal of the 
eighties [160, 161, 162, 163], which assumes the vacuum expectation value of an inflaton field as 
the source, another, perhaps more well-known, example of semiclassical gravity. 

2.1 The importance of quantum fluctuations 

For a free quantum field semiclassical gravity is fairly well understood. The theory is in some 
sense unique, since the only reasonable c-number stress-energy tensor that one may construct 
[152, 153] with the stress-energy operator is a renormalized expectation value. However the scope 
and limitations of the theory are not so well understood. It is expected that the semiclassical 
theory would break down at the Planck scale. One can conceivably assume that it would also break 
down when the fluctuations of the stress-energy operator are large [164, 165]. Calculations of the 
fluctuations of the energy density for Minkowski, Casimir and hot flat spaces as well as Einstein 
and de Sitter universes are available [165, 166, 95, 96, 7, 8, 167, 6, 94, 168, 112, 170, 171]. It is 
less clear, however, how to quantify what a large fluctuation is, and different criteria have been 
proposed [165, 172, 173, 95, 96, 174, 175]. The issue of the validity of the semiclassical gravity 
viewed in the light of quantum fluctuations is summarized in our Erice lectures [3]. One can see 
the essence of the problem by the following example inspired by Ford [164]. 

Let us assume a quantum state formed by an isolated system which consists of a superposition 
with equal amplitude of one configuration of mass M with the center of mass at Xi and another 
configuration of the same mass with the center of mass at X2 . The semiclassical theory as described 
by the semiclassical Einstein equation predicts that the center of mass of the gravitational field of 
the system is centered at (Xi + X2)/2. However, one would expect that if we send a succession 
of test particles to probe the gravitational field of the above system half of the time they would 
react to a gravitational field of mass M centered at Xi and half of the time to the field centered at 
X2. The two predictions are clearly different, note that the fluctuation in the position of the center 
of masses is of the order of {Xi — ^2)^. Although this example raises the issue of how to place 
the importance of fluctuations to the mean, a word of caution should be added to the effect that 
it should not be taken too literally. In fact, if the previous masses are macroscopic the quantum 
system decoheres very quickly [46, 50] and instead of being described by a pure quantum state 



5 



it is described by a density matrix which diagonahzes in a certain pointer basis. For observables 
associated to such a pointer basis the density matrix description is equivalent to that provided by 
a statistical ensemble. The results will differ, in any case, from the semiclassical prediction. 

In other words, one would expect that a stochastic source that describes the quantum fluctua- 
tions should enter into the semiclassical equations. A significant step in this direction was made in 
Ref. [1] where it was proposed to view the back-reaction problem in the framework of an open quan- 
tum system: the quantum fields seen as the "environment" and the gravitational field as the "sys- 
tem" . Following this proposal a systematic study of the connection between semiclassical gravity 
and open quantum systems resulted in the development of a new conceptual and technical frame- 
work where (semiclassical) Einstein-Langevin equations were derived [87, 88, 85, 89, 90, 91, 92]. 
The key technical factor to most of these results was the use of the influence functional method of 
Feynman and Vernon [12] when only the coarse-grained effect of the environment on the system is of 
interest. Note that the word semiclassical put in parentheses refers to the fact that the noise source 
in the Einstein-Langevin equation arises from the quantum field, while the background spacetime 
is classical; generally we will not carry this word since there is no confusion that the source which 
contributes to the stochastic features of this theory comes from quantum fields. 

In the language of the consistent histories formulation of quantum mechanics [52, 53, 54, 55, 
56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76] for the existence 
of a semiclassical regime for the dynamics of the system one needs two requirements: The first is 
decoherence, which guarantees that probabilities can be consistently assigned to histories describing 
the evolution of the system, and the second is that these probabilities should peak near histories 
which correspond to solutions of classical equations of motion. The effect of the environment is 
crucial, on the one hand, to provide decoherence and, on the other hand, to produce both dissipation 
and noise to the system through back-reaction, thus inducing a semiclassical stochastic dynamics 
on the system. As shown by different authors [176, 43, 44, 45, 46, 47, 48, 49, 50, 51], indeed 
over a long history predating the current revival of decoherence, stochastic semiclassical equations 
are obtained in an open quantum system after a coarse graining of the environmental degrees of 
freedom and a further coarse graining in the system variables. It is expected but has not yet been 
shown that this mechanism could also work for decoherence and classicalization of the metric field. 
Thus far, the analogy could be made formally [93] or under certain assumptions, such as adopting 
the Born-Oppenheimer approximation in quantum cosmology [177, 178]. 

An alternative axiomatic approach to the Einstein-Langevin equation without invoking the 
open system paradigm was later suggested based on the formulation of a self-consistent dynamical 
equation for a perturbative extension of semiclassical gravity able to account for the lowest order 
stress-energy fluctuations of matter fields [5]. It was shown that the same equation could be derived, 
in this general case, from the influence functional of Feynman and Vernon [6]. The field equation 
is deduced via an effective action which is computed assuming that the gravitational field is a c- 
number. The important new element in the derivation of the Einstein-Langevin equation, and of the 
stochastic gravity theory, is the physical observable that measures the stress-energy fluctuations, 
namely, the expectation value of the symmetrized bi-tensor constructed with the stress-energy 
tensor operator: the noise kernel. It is interesting to note that the Einstein-Langevin equation can 
also be understood as a useful intermediary tool to compute symmetrized two-point correlations 
of the quantum metric perturbations on the semiclassical background, independent of a suitable 
classicalization mechanism [169]. 
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3 The Einstein-Langevin equation: Axiomatic approach 



In this section we introduce stochastic semiclassical gravity, or stochastic gravity for short, in an 
axiomatic way. It is introduced as an extension of semiclassical gravity motivated by the search of 
self-consistent eqiiations which describe the back-reaction of the quantum stress-energy fluctuations 
on the gravitational field [5]. 



3.1 Semiclassical gravity 

Semiclassical gravity describes the interaction of a classical gravitational field with quantum matter 
fields. This theory can be formally derived as the leading 1/N approximation of quantum gravity 
interacting with N independent and identical free quantum fields [179, 180, 181, 182] which interact 
with gravity only. By keeping the value of NG finite, where G is Newton's gravitational constant, 
one arrives at a theory in which formally the gravitational field can be treated as a c-number field 
(i.e. quantized at tree level) while matter fields are fully quantized. The semiclassical theory may 
be summarized as follows. 

Let {M,gab) be a globally hyperbolic four-dimensional spacctime manifold Ai with metric gab 
and consider a real scalar quantum field cp of mass m propagating on that manifold; we just assume 
a scalar field for simplicity. The classical action Sm for this matter field is given by the functional 

Sm[g, ^\ = -\j d^xV^ [ff^VaC^Vbc/) + (m^ + ^r) <^2] ^ (3.1) 

where Va is the covariant derivative associated to the metric gab, ^ is a coupling parameter between 
the field and the scalar curvature of the underlying spacetime R, and g = detgab- 

The field may be quantized in the manifold using the standard canonical quantization formalism 
[32, 33, 34]. The field operator in the Heisenberg representation cp is an operator- valued distribution 
solution of the Klein-Gordon equation, the field equation derived from Eq. (3.1), 

(□ - - ^R)4> = 0. (3.2) 

We may write the field operator as ^[g; x) to indicate that it is a functional of the metric gab and a 
function of the spacetime point x. This notation will be used also for other operators and tensors. 

The classical stress-energy tensor is obtained by functional derivation of this action in the usual 
way r«\x) = {2/^)5Sm/6gab, leading to 

+^ (^g"^a - V"V^ + G"^) (/)^ (3.3) 

where □ = VaV" and Gab is the Einstein tensor. With the notation T"''^[g, 0] we explicitly indicate 
that the stress-energy tensor is a functional of the metric gab and the field 4>. 

The next step is to define a stress-energy tensor operator T"'^[g;x). Naively one would replace 
the classical field (j)[g; x) in the above functional by the quantum operator (p[g; ,t), but this procedure 
involves taking the product of two distributions at the same spacetime point. This is ill-defined 
and we need a regularization procedure. There are several regularization methods which one may 
use, one is the point-splitting or point-separation regularization method [150, 151] in which one 
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introduces a point y in a neighborhood of the point x and then uses as the regulator the vector 
tangent at the point x of the geodesic joining x and y; this method is discussed for instance in 
Refs. [96, 7, 8] and in section 5. Another well known method is dimensional rcgularization in 
which one works in arbitrary n dimensions, where n is not necessarily an integer, and then uses 
as the regulator the parameter e = n — A; this method is implicitly used in this section. The 
regularized stress-energy operator using the Weyl ordering prescription, i.e. symmetrical ordering, 
can be written as 

T%] = ^{V>[y] , vVb]} + V'^'ig] 4>'[9], (3.4) 
where V'^'^lg] is the differential operator: 

■pab = _ 1 /4) c,«^n + ^ (^R"^ - V"V'') . (3.5) 

Note that if dimensional rcgularization is used, the field operator (f)[g; x) propagates in a n- 
dimensional spacetime. Once the rcgularization prescription has been introduced a regularized 
and renormalized stress-energy operator T^^[g;x) may be defined as 

f<g[y; x) = fabig; x) + Fg[g- x)!, (3.6) 

which differs from the regularized Tab[g', x) by the identity operator times some tensor counterterms 
Fg[g;x), which depend on the regulator and are local functionals of the metric, see Ref. [6] for 
details. The field states can be chosen in such a way that for any pair of physically acceptable 
states, i.e., Hadamard states in the sense of Ref. [34], \^), and \(p) the matrix element {iplT^lcp), 
defined as the limit when the regulator takes the physical value, is finite and satisfies Wald's axioms 
[33, 152]. These counterterms can be extracted from the singular part of a Schwinger-DeWitt series 
[33, 150, 151, 183]. The choice of these counterterms is not unique but this ambiguity can be 
absorbed into the renormalized coupling constants which appear in the equations of motion for the 
gravitational field. 

The semiclassical Einstein equation for the metric gab can then be written as 

Gab[g] + ^Qab - 2(a^a6 + = 87rG(f,^[y]), (3.7) 

where {T^[g\) is the expectation value of the operator T^[g,x) after the regulator takes the phys- 
ical value in some physically acceptable state of the field on {M.,gab)- Note that both the stress 
tensor and the quantum state are functionals of the metric, hence the notation. The parameters G, 
A, a and (3 are the renormalized coupling constants, respectively, the gravitational constant, the 
cosmological constant and two dimensionless coupling constants which are zero in the classical Ein- 
stein equation. These constants must be understood as the result of "dressing" the bare constants 
which appear in the classical action before renormalization. The values of these constants must be 
determined by experiment. The left hand side of Eq. (3.7) may be derived from the gravitational 
action 

Ir- K + aCabcdC'''""' + (3R^ , (3.8) 



where Cabcd is the Weyl tensor. The tensors Aab and Bab come from the functional derivatives with 
respect to the metric of the terms quadratic in the curvature in Eq. (3.8), they are explicitly given 
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by 



^Sgab 

-2UR°^ + |v"V''ii; + \g°^^R, (3.9) 

B'^ = -^-^fd'^R' 

V-9 ogab J 

= ]^g°'^R^ - 2RR'''' + 2V''V''R - 2g"-''nR, (3.10) 

where Rabcd and Rab are the Riemann and Ricci tensors, respectively. These two tensors are, hke 
the Einstein and metric tensors, symmetric and divergenceless: V^Aab = = V^Bab- 

A sohition of semiclassical gravity consists of a spacetime {M,gab)i a quantum field operator 
which satisfies the evolution equation (3.2), and a physically acceptable state \ip[g]) for this 
field, such that Eq. (3.7) is satisfied when the expectation value of the renormalized stress-energy 
operator is evaluated in this state. 

For a free quantum field this theory is robust in the sense that it is self-consistent and fairly well 
understood. As long as the gravitational field is assumed to be described by a classical metric, the 
above semiclassical Einstein equations seems to be the only plausible dynamical equation for this 
metric: the metric couples to matter fields via the stress-energy tensor and for a given quantum 
state the only physically observable c-numbcr stress-energy tensor that one can construct is the 
above renormalized expectation value. However, lacking a full quantum gravity theory the scope 
and limits of the theory are not so well understood. It is assumed that the semiclassical theory 
should break down at Planck scales, which is when simple order of magnitude estimates suggest that 
the quantum effects of gravity should not be ignored because the energy of a quantum fluctuation 
in a Planck size region, as determined by the Heisenberg uncertainty principle, is comparable to 
the gravitational energy of that fluctuation. 

The theory is expected to break down when the fluctuations of the stress-energy operator are 
large [164]. A criterion based on the ratio of the fluctuations to the mean was proposed by Kuo 
and Ford [165] (see also work via zeta-function methods [166, 171]). This proposal was questioned 
by Phillips and Hu [95, 96, 7] because it does not contain a scale at which the theory is probed or 
how accurately the theory can be resolved. They suggested the use of a smearing scale or point- 
separation distance, for integrating over the bi-tensor quantities, equivalent to a stipulation of the 
resolution level of measurements; see also the response by Ford [172, 173]. A different criterion 
is recently suggested by Anderson et al. [174, 175] based on linear response theory. A partial 
summary of this issue can be found in our Erice Lectures [3] . 



3.2 Stochastic gravity 

The purpose of stochastic gravity is to extend the semiclassical theory to account for these fluctua- 
tions in a self-consistent way. A physical observable that describes these fluctuations to lowest order 
is the noise kernel bi-tensor, which is defined through the two point correlation of the stress-energy 



9 



operator as 

Nabcd[g\x,y) = ^{{tab[9;x),tcd[9;y)}), (3.ii) 

where the curly brackets mean anticommutator, and where 

iab[9;x) = fab[g;x) - {fab[g;x)). (3.12) 

This bi-tensor can also be written Nab,c'd' [9'-, x, y), or Nab,dd' {x, y) as we do in section 5, to emphasize 
that it is a tensor with respect to the first two indices at the point x and a tensor with respect to 
the last two indices at the point j/, but we shall not follow this notation here. The noise kernel 
is defined in terms of the unrenormalized stress-tensor operator Tab[g]x) on a given background 
metric Qab, thus a regulator is implicitly assumed on the right-hand side of Eq. (3.11). However, 
for a linear quantum field the above kernel - the expectation function of a bi-tensor - is free of 
ultraviolet divergences because the regularized Tab[9\x) differs from the renormalized T^[g\x) by 
the identity operator times some tensor counterterms, see Eq. (3.6), so that in the subtraction 
(3.12) the counterterms cancel. Consequently the ultraviolet behavior of {Tab{x)T(^{y)) is the 
same as that of {Tab{x)) {Tcd{y)) j and Tab can be replaced by the renormalized operator in 
Eq. (3.11); an alternative proof of this result is given in Ref. [7, 8]. The noise kernel should 
be thought of as a distribution function, the limit of coincidence points has meaning only in the 
sense of distributions. The bi-tensor Nabcd[9\x,y), or Nabcdix,y) for short, is real and positive 
semi-definite, as a consequence of being self-adjoint. A simple proof is given in Ref. [4]. 

Once the fluctuations of the stress-energy operator have been characterized we can perturba- 
tivcly extend the scmiclassical theory to account for such fluctuations. Thus we will assume that 
the background spacetime metric gab is a solution of the semiclassical Einstein Eqs. (3.7) and we 
will write the new metric for the extended theory as gab + hab, where we will assume that hab is a 
perturbation to the background solution. The renormalized stress-energy operator and the state of 
the quantum field may now be denoted by T^[g + h] and [^[^ -|- h]), respectively, and {T^[g + h]) 
will be the corresponding expectation value. 

Let us now introduce a Gaussian stochastic tensor field ^ab[g'i x) defined by the following corre- 
lators: 

{^ab[g;x))s = 0, {Lb[g;x)^cd[g;y))s = Nabcd[g;x,y), (3.13) 

where (. . .)s means statistical average. The symmetry and positive semi-definite property of the 
noise kernel guarantees that the stochastic field tensor (,ab[g-^ x), or ^„6(.t) for short, just introduced 
is well defined. Note that this stochastic tensor captures only partially the quantum nature of the 
fluctuations of the stress-energy operator since it assumes that cumulants of higher order are zero. 

An important property of this stochastic tensor is that it is covariantly conserved in the back- 
ground spacetime V^^abig] x) = 0. In fact, as a consequence of the conservation of T^[g] one can sec 
that V^Nabcd{x,y) = 0. Taking the divergence in Eq. (3.13) one can then show that {V^^ab)s = 
and {Vx^ab{x)^cd{y)) s = so that V^^ab is deterministic and represents with certainty the zero 
vector field in Ai . 

For a conformal field, i.e., a field whose classical action is conformally invariant, is traccless: 
g"'^^ab[g'-i x) = 0; so that, for a conformal matter field the stochastic source gives no correction to the 
trace anomaly. In fact, from the trace anomaly result which states that g"'''T^[g] is, in this case, a 
local c-number functional of gab times the identity operator, we have that g"'^{x)Nabcd[g] x,y) = 0. 
It then follows from Eq. (3.13) that {g"'^£,ab)s = and {g"''^{x)£,ab{x)S,cd{y))s = 0; an alternative 
proof based on the point-separation method is given in Ref. [7, 8], see also section 5. 
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All these properties make it quite natural to incorporate into tlic Einstein equations the stress- 
energy fluctuations by using the stochastic tensor ^ab[9'i x) as the source of the metric perturbations. 
Thus we will write the following equation. 

Gab[g+h]+Ai9ab+hab) - 2{aAab + l3Bab)[g+h]=87rG({f^^[g+h])+U[9i) ■ (3.14) 

This equation is in the form of a (semiclassical) Einstein-Langevin equation, it is a dynamical 
equation for the metric perturbation hab to linear order. It describes the back-reaction of the 
metric to the quantum fluctuations of the stress-energy tensor of matter fields, and gives a first 
order extension to semiclassical gravity as described by the semiclassical Einstein equation (3.7). 

Note that we refer to the Einstein-Langevin equation as a first order extension to semiclassical 
Einstein equation of semiclassical gravity and the lowest level representation of stochastic gravity. 
However, stochastic gravity has a much broader meaning, it refers to the range of theories based on 
second and higher order correlation functions. Noise can be defined in effectively open systems (e.g. 
correlation noise [29] in the Schwinger-Dyson equation hierarchy) to some degree but one should 
not expect the Langevin form to prevail. In this sense we say stochastic gravity is the intermediate 
theory between semiclassical gravity (a mean field theory based on the expectation values of the 
energy momentum tensor of quantum fields) and quantum gravity (the full hierarchy of correlation 
functions retaining complete quantum coherence [2, 133]. 

The renormalization of the operator Tab[g + h] is carried out exactly as in the previous case, 
now in the perturbed metric gab + hab- Note that the stochastic source $,ab[g',x) is not dynamical, 
it is independent of hab since it describes the fluctuations of the stress tensor on the semiclassical 
background gab- 

An important property of the Einstein-Langevin equation is that it is gauge invariant under the 
change of hab by /i^^ = hab + ^aCb + ^bCa, where is a stochastic vector field on the background 
manifold M- Note that a tensor such as Rab[g + h], transforms as Rab[g + h'] = Rab[g + h] +Ci^Rab[g] 
to linear order in the perturbations, where is the Lie derivative with respect to C"- Now, let 
us write the source tensors in Eqs. (3.14) and (3.7) to the left-hand sides of these equations. If we 
substitute h by h' in this new version of Eq. (3.14), we get the same expression, with h instead of h' , 
plus the Lie derivative of the combination of tensors which appear on the left-hand side of the new 
Eq. (3.7). This last combination vanishes when Eq. (3.7) is satisfied, i.e., when the background 
metric gab is a solution of semiclassical gravity. 

A solution of Eq. (3.14) can be formally written as hab[^]- This solution is characterized by 
the whole family of its correlation functions. From the statistical average of this equation we 
have that gab + {hab)s must be a solution of the semiclassical Einstein equation linearized around 
the background gab', this solution has been proposed as a test for the validity of the semiclassical 
approximation [174, 175]. The fluctuation of the metric around this average are described by the 
moments of the stochastic field h^abiH ~ habiCi — {hab)s- Thus the solutions of the Einstein-Langevin 
equation will provide the two-point metric correlation functions {h^jj{x)h^j^{y))s- 

We see that whereas the semiclassical theory depends on the expectation value of the point- 
defined value of the stress-energy operator, the stochastic theory carries information also on the two 
point correlation of the stress-energy operator. We should also emphasize that, even if the metric 
fluctuations appears classical and stochastic, their origin is quantum not only because they are 
induced by the fluctuations of quantum matter, but also because they are the suitably coarse-grained 
variables left over from the quantum gravity fluctuations after some mechanism for decoherence 
and classicalization of the metric field [176, 184, 61, 62, 63, 185]. One may, in fact, derive the 
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stochastic semiclassical theory from a full quantum theory. This was done via the world-line 
influence functional method for a moving charged particle in an electromagnetic field in quantum 
electrodynamics [186]. From another viewpoint, quite independent of whether a classicalization 
mechanism is mandatory or implementable, the Einstein-Langevin equation proves to be a useful 
tool to compute the symmetrized two point correlations of the quantum metric perturbations [169]. 
This is illustrated in the linear toy model discussed in Ref. [4], which has features of some quantum 
Brownian models [187, 188, 189]. 

4 The Einstein-Langevin equation: Functional approach 

The Einstein-Langevin equation (3.14) may also be derived by a method based on functional tech- 
niques [6]. Here we will summarize these techniques starting with semiclassical gravity. 

In semiclassical gravity functional methods were used to study the back-reaction of quantum 
fields in cosmological models [154, 155, 156]. The primary advantage of the effective action approach 
is, in addition to the well-known fact that it is easy to introduce perturbation schemes like loop 
expansion and nPI formalisms, that it yields a fully self-consistent solution. For a general discussion 
in the semiclassical context of these two approaches, equation of motion versus effective action, see, 
e.g., the work of Hu and Parker (1978) versus Hartle and Hu (1979) in [36, 37, 38, 39, 40, 41, 42]. 
See also comments in Sec. 5.6 on the black hole backreaction problem comparing the approach by 
York et al. [190, 191, 192] versus that of Sinha, Raval and Hu [123]. 

The well known in-out effective action method treated in textbooks, however, led to equations of 
motion which were not real because they were tailored to compute transition elements of quantum 
operators rather than expectation values. The correct technique to use for the backreaction problem 
is the Schwinger-Keldysh [14, 15, 16, 17, 18, 19, 20] closed-time-path (CTP) or 'in-in' effective 
action. These techniques were adapted to the gravitational context [21, 22, 23, 24, 25, 84] and 
applied to different problems in cosmology. One could deduce the semiclassical Einstein equation 
from the CTP effective action for the gravitational field (at tree level) with quantum matter fields. 

Furthermore, in this case the CTP functional formalism turns out to be related [18, 87, 89, 
193, 194, 121, 122, 195, 196, 6, 93] to the influence functional formalism of Feynman and Vernon 
[12] since the full quantum system may be understood as consisting of a distinguished subsystem 
(the "system" of interest) interacting with the remaining degrees of freedom (the environment). 
Integrating out the environment variables in a CTP path integral yields the influence functional, 
from which one can define an effective action for the dynamics of the system [87, 85, 197, 194]. This 
approach to semiclassical gravity is motivated by the observation [1] that in some open quantum 
systems classicalization and decoherence [43, 44, 45, 46, 47, 48, 49, 50, 51] on the system may be 
brought about by interaction with an environment, the environment being in this case the matter 
fields and some "high-momentum" gravitational modes [77, 78, 79, 80, 81, 82, 83, 185]. Unfortu- 
nately, since the form of a complete quantum theory of gravity interacting with matter is unknown, 
we do not know what these "high-momentum" gravitational modes are. Such a fundamental quan- 
tum theory might not even be a field theory, in which case the metric and scalar fields would not 
be fundamental objects [2]. Thus, in this case, we cannot attempt to evaluate the influence action 
of Feynman and Vernon starting from the fundamental quantum theory and performing the path 
integrations in the environment variables. Instead, we introduce the influence action for an effec- 
tive quantum flcld theory of gravity and matter [198, 199, 200, 201, 202, 177, 178], in which such 
"high-momentum" gravitational modes are assumed to have already been "integrated out." 



12 



4.1 Influence action for semiclassical gravity 

Let us formulate semiclassical gravity in this functional framework. Adopting the usual procedure 

of effective field theories [203, 204, 198, 199, 200, 201, 205], one has to take the effective action for 
the metric and the scalar field of the most general local form compatible with general covariance: 
S[g,(l)] = Sg[g] + Sm[g,<f)\ + where Sg[g] and Sm[g,<p] are given by Eqs. (3.8) and (3.1), 
respectively, and the dots stand for terms of order higher than two in the curvature and in the 
number of derivatives of the scalar field. Here, we shall neglect the higher order terms as well 
as self-interaction terms for the scalar field. The second order terms are necessary to renormalize 
one-loop ultraviolet divergences of the scalar field stress-energy tensor, as we have already seen. 
Since M. is a globally hyperbolic manifold, we can foliate it by a family of t = constant Cauchy 
hypersurfaces S^, and we will indicate the initial and final times by tj and tf, respectively. 

The influence functional corresponding to the action (3.1) describing a scalar field in a spacetime 
(coupled to a metric field) may be introduced as a functional of two copies of the metric, denoted 
by g^f^ and g^j^, which coincide at some final time t = tf. Let us assume that, in the quantum 
effective theory, the state of the full system (the scalar and the metric fields) in the Schrodinger 
picture at the initial time t = ti can be described by a density operator which can be written as 
the tensor product of two operators on the Hilbert spaces of the metric and of the scalar field. Let 
Pi{ti) = Pi [0+(ii), 0_(ti)] be the matrix element of the density operator p^{ti) describing the initial 
state of the scalar field. The Feynman- Vernon influence functional is defined as the following path 
integral over the two copies of the scalar field: 

J^jF[g^] = J V<t>+ V<p_ p,{t,)6[Mtf)-Mtf)] e^^-I^+'-^+l-^-t^-'^-l). (4.1) 

Alternatively, the above double path integral can be rewritten as a closed time path (CTP) integral, 
namely, as a single path integral in a complex time contour with two different time branches, one 
going forward in time from ti to tf, and the other going backward in time from tf to ti (in practice 
one usually takes ti — > — oo). Prom this influence functional, the influence action Swlg'^ , g~], or 
•SiFb^] for short, defined by 

J^iF[5±] = e^^i^[^*], (4.2) 

carries all the information about the environment (the matter fields) relevant to the system (the 
gravitational field). Then we can define the CTP effective action for the gravitational field, >S'eff [5^], 
as 

S,s[g^] = Sg[g+] - Sg[g-] + SMg^]. (4.3) 

This is the effective action for the classical gravitational field in the CTP formalism. However, 
since the gravitational field is treated only at the tree level, this is also the effective classical action 
from which the classical equations of motion can be derived. 

Expression (4.1) contains ultraviolet divergences and must be regularized. We shall assume 
that dimensional regularization can be applied, that is, it makes sense to dimensionally continue all 
the quantities that appear in Eq. (4.1). For this we need to work with the n-dimensional actions 
corresponding to Sm in (4.1) and Sg in (3.8). For example, the parameters G, A a and (3 of Eq. (3.8) 
are the bare parameters Gb, A^, ob and Pb, and in Sg[g], instead of the square of the Weyl tensor 
in Eq. (3.8), one must use {2/3){RabcdR"'^"'' — RabR""^) which by the Gauss-Bonnet theorem leads 
to the same equations of motion as the action (3.8) when n = 4. The form of Sg in n dimensions 
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is suggested by the Schwinger-DeWitt analysis of the ultraviolet divergences in the matter stress- 
energy tensor using dimensional regularization. One can then write the Feynman- Vernon effective 
action [5^] in Eq. (4.3) in a form suitable for dimensional regularization. Since both Sm and Sg 
contain second order derivatives of the metric, one should also add some boundary terms [31, 85]. 
The effect of these terms is to cancel out the boundary terms which appear when taking variations 
of ScQ[g^] keeping the value of g'^f^ and g~f^ fixed at E^. and Alternatively, in order to obtain 
the equations of motion for the metric in the semiclassical regime, we can work with the action 
terms without boundary terms and neglect all boundary terms when taking variations with respect 
to g"^^. Prom now on, all the functional derivatives with respect to the metric will be understood 
in this sense. 

The semiclassical Einstein equation (3.7) can now be derived. Using the definition of the stress- 
energy tensor r"^(x) = {2/ y/^)5Sm/ ^9ah and the definition of the influence functional, Eqs. (4.1) 
and (4.2), we see that 



(4.4) 

=9 

where the expectation value is taken in the n-dimensional spacetime generalization of the state 
described by p^{ti). Therefore, differentiating S'egf^^] in Eq. (4.3) with respect to g'^jj, and then 
setting = g~j^ = gab, we get the semiclassical Einstein equation in n dimensions. This equation is 
then renormalized by absorbing the divergences in the regularized {T°'^[g]) into the bare parameters. 
Taking the limit n — 4 we obtain the physical semiclassical Einstein equation (3.7). 

4.2 Influence action for stochastic gravity 

In the spirit of the previous derivation of the Einstein-Langevin equation, we now seek a dynamical 
equation for a linear perturbation hab to the semiclassical metric gab, solution of Eq. (3.7). Strictly 
speaking if we use dimensional regularization we must consider the n-dimensional version of that 
equation. From the results just described, if such an equation were simply a linearized semiclassical 
Einstein equation, it could be obtained from an expansion of the effective action Ses[g + h^]. In 
particular, since, from Eq. (4.4), we have that 



V'-det((7+/i)(x) Sh+f,{. 



(4.5) 

h^=h 



the expansion of {f"'^ [g+h] ) to linear order in hab can be obtained from an expansion of the influence 
action SiF[g + h^] up to second order in h^. 

To perform the expansion of the influence action, we have to compute the flrst and second order 
functional derivatives of SiF[g + h^] and then set h'^b~^ab~^0'b- If we do so using the path integral 
representation (4.1), wc can interpret these derivatives as expectation values of operators. The 
relevant second order derivatives are 



6^SiF[g + h^ 



4 S^Swig^] 



= -H^'"^[g-x,y)-K'^'^%-x,y) 

h^=h 



V-9{x)V-9{y) 5Kix)6h^ly) 



= -Hf^[g; X, y)-iN'''"^%- x, y\ (4.6) 

/i±=/i 
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where 

N'^^'^\g-x,y)^\{{i'^\g-x)J'%-y)} 
Hf<^%; X, y) = Im {T*(f-%; x)f-%- y)) ), 
Hf<^[g-x,y) ^ ([r«^b;x), f'^%-.y)^. 

K-^-%-X,y) = / ^Srr^g + h^ 

with i""^ defined in Eq. (3.12), [ , ] denotes the commutator and { , } the anti-commutator. Here we 
use a Weyl ordering prescription for the operators. The symbol T* denotes the following ordered 
operations: First, time order the field operators (f) and then apply the derivative operators which 
appear in each term of the product T""^ {x)T^'^ {y) , where T"'' is the functional (3.3). This T* "time 
ordering" arises because we have path integrals containing products of derivatives of the field, which 
can be expressed as derivatives of the path integrals which do not contain such derivatives. Notice, 
from their definitions, that all the kernels which appear in expressions (4.6) are real and also H^'^ 
is free of ultraviolet divergences in the limit n ^ 4. 

From (4.4) and (4.6), since SiY[g-,g\ = and SiY[g~ig'^] = — S^^ig'^ , g~\-, we can write the 
expansion for the influence action Sip [g + around a background metric gab in terms of the 
previous kernels. Taking into account that these kernels satisfy the symmetry relations 

Ht'^'^ix, y) = Hf'^\y, x),Hf'''ix, y) = -Hf'^\y, x),K''''^\x, y)=K^''\y, x), (4.7) 

and introducing the new kernel 

H'^'"^{x, y) = Hf'^ix, y) + Hf'^{x, y), (4.8) 

the expansion of ^ip can be finally written as 

SMg+h^] = \jd!'x,p^){f-%-x))[hab{x)\ 
j d^xd'^y ^J-g{x)sJ-g{y) 

X [Kb{x)\ (H'^'^[g-x,y)+K'^'"^%-x,y)) {hM) 

d^xd^y sj -g{x)^J -g{y) 
X [h^{x)] X, y) [hM + 0(/i'), (4.9) 

where we have used the notation 

\Kb\ ^ h+,-h-„ {Kb} ^ h++h-,. (4.10) 

From Eqs. (4.9) and (4.5) it is clear that the imaginary part of the influence action does not 
contribute to the perturbed semiclassical Einstein equation (the expectation value of the stress- 
energy tensor is real), however, as it depends on the noise kernel, it contains information on the 
fluctuations of the operator T^-^lg]. 



15 



We are now in a position to carry out the derivation of the semiclassical Einstein-Langevin 
equation. The procedure is well known [87, 85, 89, 206, 207, 208, 209]: it consists of deriving a 
new "stochastic" effective action from the observation that the effect of the imaginary part of the 
influence action (4.9) on the corresponding influence functional is equivalent to the averaged effect 
of the stochastic source coupled linearly to the perturbations K^. This observation follows from 
the identity first invoked by Feynman and Vernon for such purpose: 

exp (^-ljd^xd%^^)^^) [hab{x)\ N'^'''\x,y) [KM]) 

= jviVm exv{^i^jd^x,p^)e\x) [hab{x)]) , (4.11) 

where is the probability distribution functional of a Gaussian stochastic tensor characterized 
by the correlators (3.13) with JSl"-^'^'^ given by Eq. (3.11), and where the path integration measure 
is assumed to be a scalar under diffeomorphisms of {M.,gah)- The above identity follows from the 
identification of the right-hand side of (4.11) with the characteristic functional for the stochastic 
field The probability distribution functional for is explicitly given by 



V[i] = dei{2TTN)-^/\-x.^ 



(4.12) 



2. 

We may now introduce the stochastic effective action as 

Slw[g + h^,^] = Sg[g + h+] - Sg[g + + S^g + h^,^], (4.13) 
where the "stochastic" influence action is defined as 

SfAg + h^,^] = ^Sw[g+h^] +^ J d^x ^J-g{x)e\x) [hab{x)] + 0{h^). (4.14) 

Note that, in fact, the influence functional can now be written as a statistical average over 
^w[g + h^] = (exp {iSfplg + ,i\))s ■ The stochastic equation of motion for hab reads 



5Sls[g + h^,£] 



= 0, (4.15) 

h^=h 



which is the Einstein-Langevin equation (3.14); notice that only the real part of Sjf contributes 

to the expectation value (4.5). To be precise we get first the regularized n-dimensional equations 
with the bare parameters, and where instead of the tensor we get (2/3)L'"* where 



1 5 



-g ^9. 



ab 



-g (RcdefR'"'^ - RcdR"") 



= \ ^^'^ {RcdefR"^^^ - RcdR^ + or) - 2R'''^^R\de 

-2R°^^'^Rcd + AR°^Rc^ - 3ni?"^ + VV^R. (4.16) 

Of course, when n = 4 these tensors are related, A"^ = (2/3) L'"^ After that we renormalize and 
take the limit n — > 4 to obtain the Einstein-Langevin equations in the physical spacetime. 
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4.3 Explicit form of the Einstein-Langevin equation 



We can write the Einstein-Langevin equation in a more explicit form by working out the expansion 
of {T°'^[g+h]) up to Hnear order in the perturbation hab- Prom Eq. (4.5), we see that this expansion 
can be easily obtained from (4.9). The result is 



-^d-y ^f^)Hf"'[9; X, y)Ka{y) + 0{h^). 



(4.17) 



Here we use a subscript n on a given tensor to indicate that we are explicitly working in n- 
dimensions, as we use dimensional regularization, and we also use the superindex (i) to generally 
indicate that the tensor is the first order correction, linear in hab, in a perturbative expansion 
around the background gab- 

Using the Klein-Gordon equation (3.2), and expressions (3.3) for the stress-energy tensor and 
the corresponding operator, wc can write 



Ti^'^'ig, h]= - g^%a - S'X - S'Xd T^M + -^"'b, h] ^l[g], 



(4.18) 



where J-'°'''[g; h] is the differential operator 



h' 



ab 



(4.19) 



It is Tindcrstood that indices are raised with the background inverse metric g""^ and that all the 
covariant derivatives arc associated to the metric ^'a^. 

Substituting (4.17) into the 77,-dimensional version of the Einstein-Langevin Eq. (3.14), taking 
into account that gab satisfies the semiclassical Einstein equation (3.7), and substituting expression 
(4.18) we can write the Einstein-Langevin equation in dimensional regularization as 



1 



1 



ab /~<cd I 



^acub 



bci,a 



,ab 



„abuc 



SttGb 

- ^ (^D^nab _ ^gabj^cd^^^ ^ jjac^b ^ jjbc^a 

- K\^l[9; x)) + \j X, y)hM 



-(n-4)pab 
Sn ' 



(4.20) 



where the tensors G"*, D""^ and B"'^ are computed from the semiclassical metric gab, and where 
we have omitted the functional dependence on gab and hab in G*^^"^, D*^'"^, j^W"-b ^^^^ j:ab 
simplify the notation. The parameter is a mass scale which relates the dimensions of the physical 
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field (p with the dimensions of the corresponding field in n-dimensions, cpn = Notice 
that, in Eq. (4.20), all the ultraviolet divergences in the limit n ^ 4, which must be removed 
by renormalization of the coupling constants, are in and the symmetric part H^^'^'^(x,y) 

of the kernel H^^'^'^(x,y), whereas the kernels N^'^'^{x,y) and H^^^'^{x,y) are free of ultraviolet 
divergences. If we introduce the bi-tensor F^'"^'^[g;x,y) defined by 

Ff^^; X, y) ^ (ttig; x) iP-[g; y)) (4.21) 

where t"'^ is defined by Eq. (3.12), then the kernels N and Ha can be written as 

K''^%x,y) = ReF:'^'^%x,y), H^t%;x,y) = ImFf '^'^[y; x, y), (4.22) 

where we have used that 2{t''''{x) i'^{y)) = ({^^(x), t^'^(y)}) + ([r*(x), t'^'^iy)]), and the fact that the 
first term on the right hand side of this identity is real, whereas the second one is pure imaginary. 
Once we perform the renormalization procedure in Eq. (4.20), setting n = 4 will yield the physical 
Einstein-Langevin equation. Due to the presence of the kernel H^'^'^ix^y), this equation will be 
usually non-local in the metric perturbation. In section 6 we will carry out an explicit evaluation 
of the physical Einstein-Langevin equation which will illustrate the procedure. 



4.3.1 The kernels for the vacuum state 

When the expectation values in the Einstein-Langevin equation are taken in a vacuum state |0), 
such as, for instance, an "in" vacuum, we can be more explicit, since we can write the expectation 
values in terms of the Wightman and Feynman functions, defined as 

G+[g; X, y) = (O^^^b; x)0nb; y)|0), zG^^b; x, y) = (O[r(0„[5; x)^^; y)) |0). (4.23) 

These expressions for the kernels in the Einstein-Langevin equation will be very useful for explicit 
calculations. To simplify the notation, we omit the functional dependence on the semiclassical 
metric gab, which will be understood in all the expressions below. 

Prom Eqs. (4.22), we see that the kernels Nf'"^{x, y) and Hf'"^{x, y) are the real and imaginary 
parts, respectively, of the bi-tensor F^'^'^(x,y). Prom the expression (3.4) we see that the stress- 
energy operator can be written as a sum of terms of the form ^Ax4>n{x), 6^0„(x)|, where Ax 

and Bx are some differential operators. It then follows that we can express the bi-tensor F^'^'^{x,y) 
in terms of the Wightman function as 

F^'^^^x^y) = V«V^G+(x,y)VtV^G+(x,y) + VSV^G+(x, 2/)V^V^G+(x, y) 
+ 2Pf(V^G+(x,y)V^G+(x,y)) 

+2Vl\y%GUx,y)VlG+{x,y)) + 2VfVl\G+\x,y)), (4.24) 

where V"^ is the differential operator (3.5). Prom this expression and the relations (4.22), we get 
expressions for the kernels Nn and H^^ in terms of the Wightman function G^(x,y). 
Similarly the kernel H^^'^'^{x, y), can be written in terms of the Peynman function as 



H^^'\x,y) = -Im VSV^G,Jx,y)V^V^ajx,y) 



b\-7d. 
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+V«V^a.„(x,y)V^V^a.„(x,y) 

- 9'^\x)VlVlG,^{x, y) V/V^G,Jx, y) 

-g''\y)V%VlG,^{x, y) V^V/G,„(x, y) 

+/Cf(2V^,G^„(.T,y)V^G^Jx,y) 

^""(x) V^O.Jx, y) V/G,„(x, y)) +2/Cf /C^-^tC^Cx, y))] , (4.25) 

where K,'^ is the differential operator 

/Cf ^ i (f\x)n, - V^V^ + G-\x)) - \ m'g-\x). (4.26) 

Note that, in the vacuum state |0), the term (^^(x)) in equation (4.20) can also be written as 
(<^n(a^)) = iGi.„(x,,T) = G+(x,x). 

Finally, the causality of the Einstein-Langevin equation (4.20) can be explicitly seen as follows. 
The non-local terms in that equation are due to the kernel H{x,y) which is defined in Eq. (4.8) as 
the sum of Hs{x, y) and Ha{x, y). Now, when the points x and y are spacelike separated, (t>n{x) and 
4>n{y) commute and, thus, G^(x, y) = iGp„(x, y) = (1/2)(0| {(/!)„(x), (?!>„(y)} |0), which is real. Hence, 
from the above expressions, we have that H^'^'^'^{x,y) = H^^'^'^ (x , y) = 0, and thus Hj^"\x.y) = 0. 
This fact is expected since, from the causality of the expectation value of the stress-energy operator 
[152], we know that the non-local dependence on the metric perturbation in the Einstein-Langevin 
equation, see Eq. (3.14), must be causal. See Ref. [4] for an alternative proof of the causal nature 
of the Einstein-Langevin equation. 



5 Noise Kernel and Point- Separation 

In this section we explore further the properties of the noise kernel and the stress energy bi-tensor. 
Similar to what was done for the stress energy tensor it is desirable to relate the noise kernel defined 
at separated points to the Green function of a quantum field. We pointed out earlier [2] that field 
quantities defined at two separated points may possess important information which could be the 
starting point for probes into possible extended structures of spacetime. Of more practical concern 
is how one can define a finite quantity at one point or in some small region around it from the noise 
kernel defined at two separated points. When we refer to, say, the fiuctuations of energy density 
in ordinary (point-wise) quantum field theory, we are in actuality asking such a question. This is 
essential for addressing fundamental issues like a) the validity of semiclassical gravity [165] - e.g., 
whether the fluctuations to mean ratio is a correct criterion [95, 96, 172, 173, 174, 175]; b) Whether 
the fluctuations in the vacuum energy density which drives some models of inflationary cosmology 
violates the positive energy condition; c) Physical effects of black hole horizon fluctuations and 
Hawking radiation backreaction - to begin with, is the fluctuations flnite or infinite? d) General 
relativity as a low energy effective theory in the geometro-hydrodynamic limit towards a kinetic 
theory approach to quantum gravity [131, 2, 133]. 
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Thus, for comparison with ordinary phenomena at low energy we need to find a reasonable 
prescription for obtaining a finite quantity of the noise kernel in the limit of ordinary (point- 
defined) quantum field theory. Regularization schemes used in obtaining a finite expression for 
the stress energy tensor have been applied to the noise kernel^. This includes the simple normal 
ordering [165, 210] and smeared field operator [96] methods applied to the Minkowski and Casimir 
spaces, zeta-function [211, 212, 213] for spacetimes with an Euclidean section, applied to the Casimir 
effect [171] and the Einstein Universe [166], or the covariant point-separation methods applied to 
the Minkowski [96], hot fiat space and the Schwarzschild spacetime [8]. There are differences and 
deliberations on whether it is meaningful to seek a point- wise expression for the noise kernel, and if 
so what is the correct way to proceed - e.g., regularization by a subtraction scheme or by integrating 
over a test-field. Intuitively the smear field method [96] may better preserve the integrity of the noise 
kernel as it provides a sampling of the two point function rather than using a subtraction scheme 
which alters its innate properties by forcing a nonlocal quantity into a local one. More investigation 
is needed to clarify these points, which bear on important issues like the validity of semiclassical 
gravity. We shall set a more modest goal here, to derive a general expression for the noise kernel for 
quantum fields in an arbitrary curved spacetime in terms of Green functions and leave the discussion 
of point-wise limit to a later date. For this purpose the covariant point-separation method which 
highlights the bi-tensor features, when used not as a regularization scheme, is perhaps closest to 
the spirit of stochastic gravity. 

The task of finding a general expression of the noise-kernel for quantum fields in curved space- 
times was carried out by Phillips and Hu in two papers using the "modified" point separation 
scheme [149, 214, 153]. Their first paper [7] begins with a discussion of the procedures for dealing 
with the quantum stress tensor bi-operator at two separated points, and ends with a general expres- 
sion of the noise kernel defined at separated points expressed as products of covariant derivatives 
up to the fourth order of the quantum field's Green function. (The stress tensor involves up to two 
covariant derivatives.) This result holds for x ^ y without recourse to renormalization of the Green 
function, showing that Nabc'd'{x,y) is always finite for x ^ y (and off the light cone for massless 
theories). In particular, for a massless conformally coupled free scalar field on a four dimensional 
manifold they computed the trace of the noise kernel at both points and found this double trace 
vanishes identically. This implies that there is no stochastic correction to the trace anomaly for 
massless conformal fields, in agreement with results arrived at in Refs. [87, 89, 6] (see also section 
3). In their second paper [8] a Gaussian approximation for the Green function (which is what limits 
the accuracy of the results) is used to derive finite expressions for two specific classes of spacetimes, 
ultrastatic spacetimes, such as the hot flat space, and the conformally- ultrastatic spacetimes, such 
as the Schwarzschild spacetime. Again, the validity of these results may depend on how we view 
the relevance and meaning of regularization. We will only report the result of their first paper here. 

^It is well-known that several regularization methods can work equally well for the removal of ultraviolet diver- 
gences in the stress energy tensor of quantum fields in curved spacetime. Their mutual relations are known, and 
discrepancies explained. This formal structure of regularization schemes for quantum fields in curved spacetime 
should remain intact when applied to the regularization of the noise kernel in general curved spacetimes; it is the 
meaning and relevance of regularization of the noise kernel which is more of a concern (see comments below) . Specific 
considerations will of course enter for each method. But for the methods employed so far, such as zeta-function, 
point separation, dimensional, smeared-field, applied to simple cases (Casimir, Einstein, thermal fields) there is no 
new inconsistency or discrepancy. 
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5.1 Point Separation 



The point separation scheme introduced in the 60's by DeWitt [215] was brought to more popular 

use in the 70's in the context of quantum field theory in curved spacetimes [134, 150, 151] as a 
means for obtaining a finite quantum stress tensor. Since the stress-energy tensor is buih from 
the product of a pair of field operators evaluated at a single point, it is not well-defined. In this 
scheme, one introduces an artificial separation of the single point x to a pair of closely separated 
points X and x'. The problematic terms involving field products such as becomes (f){x)(f){x'), 

whose expectation value is well defined. If one is interested in the low energy behavior captured 
by the point-defined quantum field theory - as the effort in the 70's was directed - one takes the 
coincidence limit. Once the divergences present are identified, they may be removed (regularization) 
or moved (by renormalizing the coupling constants), to produce a well-defined, finite stress tensor 
at a single point. 

Thus the first order of business is the construction of the stress tensor and then derive the 
symmetric stress-energy tensor two point function, the noise kernel, in terms of the Wightman 
Green function. In this section we will use the traditional notation for index tensors in the point- 
separation context. 



5.1.1 n-tensors and end-point expansions 



An object like the Green function G{x, y) is an example of a bi-scalar. it transforms as scalar 
at both points x and y. We can also define a bi-tensor r^^...^^ (x, y): upon a coordinate 
transformation, this transforms as a rank n tensor at x and a rank m tensor at y. We will extend 
this up to a quad-tensor T„ b' ■■■b' c"---c" d"'---d"' which has support at four points x,y,x',y', 
transforming as rank ni,n2, 713,714 tensors at each of the four points. This also sets the notation 
we will use: unprimed indices referring to the tangent space constructed above x, single primed 
indices to y, double primed to x' and triple primed to y'. For each point, there is the covariant 
derivative Va at that point. Covariant derivatives at different points commute and the covariant 
derivative at, say, point x', does not act on a bi-tensor defined at, say, x and y: 



= Tab';d';c and Tab'-c" = 0- 



(5.1) 



To simplify notation, henceforth we will eliminate the semicolons after the first one for multiple 
covariant derivatives at multiple points. 

Having objects defined at different points, the coincident limit is defined as evaluation "on the 
diagonal" , in the sense of the spacetime support of the function or tensor, and the usual shorthand 
[G{x, y)] = G{x, x) is used. This extends to 7T,-tensors as 



h' ■■■h' r"---r" d"'---d' 



■*n4 



'ai---ani b\—bn2 Cl---Cn^ di--d„^ ■ 



(5.2) 



i.e., this becomes a rank (rzi + n2 + + n^) tensor at x. The multi- variable chain rule relates 
covariant derivatives acting at different points, when we are interested in the coincident limit: 



T, 



+ 



'^ai—amb[--b'n;c' 



(5.3) 



This result is referred to as Synge's theorem in this context; we follow Fulling's [33] discussion. 
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The bi-tensor of parallel transport ga^ is defined such that when it acts on a vector Vb' at y, it 
parallel transports the vector along the geodesies connecting x and y. This allows us to add vectors 
and tensors defined at different points. We cannot directly add a vector Va at x and vector Wa' at 
y. But by using gj' , we can construct the sum + gj' wy. We will also need the obvious property 
a - o ^ 

i/a — ya ■ 

The main bi-scalar we need is the world function a{x,y). This is defined as a half of the square 
of the geodesic distance between the points x and y. It satisfies the equation 

a=^a'Pa,, (5.4) 

Often in the literature, a covariant derivative is implied when the world function appears with 
indices: a"' = a'", i.e. taking the covariant derivative at x, while a"" means the covariant derivative 
at y. This is done since the vector —a"- is the tangent vector to the geodesic with length equal 
the distance between x and y. As (t" records information about distance and direction for the two 
points this makes it ideal for constructing a series expansion of a bi-scalar. The end point expansion 
of a bi-scalar S{x, y) is of the form 

S{x, y) = + a^'AW + aVMg + aW-^g, + aWV^^g,, + • • • (5.5) 

where, following our convention, the expansion tensors A^a\-an with unprimed indices have support 
at X and hence the name end point expansion. Only the symmetric part of these tensors contribute 
to the expansion. For the purposes of multiplying series expansions it is convenient to separate the 
distance dependence from the direction dependence. This is done by introducing the unit vector 
pa _ ^aj ,j2a. Then the series expansion can be written 

S{x, y) = ^(°) + + + ^§^(3) ^ ^2^(4) ^ . . . ^^ g^ 

The expansion scalars are related, via 

^pPi . . -pP", to the expansion tensors. 
The last object we need is the VanVleck-Morette determinant D{x,y), defined as D(x,y) = 
— det {—o'-ab')- The related bi-scalar 

satisfies the equation 

(4 - a./) - 2AV2 .^a'P = (5.8) 



with the boundary condition 



Ai/2' 



1. 

Further details on these objects and discussions of the definitions and properties are contained 
in [150, 151] and [216]. There it is shown how the defining equations for a and ^^"^ are used to 
determine the coincident limit expression for the various covariant derivatives of the world function 
([(T;a], [c-^ab], ctc.) and how the defining differential equation for A^/^ can be used to determine 
the series expansion of A^/^. We show how the expansion tensors A^]..a„ are determined in terms 
of the coincident limits of covariant derivatives of the bi-scalar S{x,y). (Ref. [216] details how 
point separation can be implemented on the computer to provide easy access to a wider range of 
applications involving higher derivatives of the curvature tensors.) 
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5.2 Stress Energy Bi- Tensor Operator and Noise Kernel 

Even though we believe that the point-separated results are more basic in the sense that it reflects 
a deeper structure of the quantum theory of spacetime, we will nevertheless start with quantities 
defined at one point because they are what enter in conventional quantum field theory. We will use 
point separation to introduce the bi-quantities. The key issue here is thus the distinction between 
point-defined {pt) and point-separated (bi) quantities. 

For a free classical scalar field <p with the action Sm[g,4>] defined in Eq. (3.1), the classical 
stress-energy tensor is 

Tab = (1-20 <P;a<P;b+(2^- \^ <P;p ^P'" 9ab + 2^ (</>;/ - <j>-ab gab) 

+02 e (Uab - ^Rgab^ - gab, (5.9) 

which is equivalent to the tensor of Eq. (3.3) but written in a slightly different form for convenience. 
When we make the transition to quantum field theory, we promote the field (j){x) to a field operator 
(f){x). The fundamental problem of defining a quantum operator for the stress tensor is immediately 
visible: the field operator appears quadratically. Since ^(a;) is an operator-valued distribution, 
products at a single point axe not well-defined. But if the product is point separated ((^^(x) 
(f>{x)(f>{x')), they are finite and well-defined. 

Let us first seek a point-separated extension of these classical quantities and then consider 
the quantum field operators. Point separation is symmetrically extended to products of covariant 
derivatives of the field according to 

{ct>;a) {ct>;b) ^ \ {ga^'^p'^b + gb^'^a^p) <t>{x)<t>{x'), 
<t>{4>;ab) ^ ^(V„Vfe + 5/VVp.V,.)</.(x)</.(x'). 

The bi- vector of parallel displacement ga"'' {x,x') is included so that we may have objects that are 
rank 2 tensors at x and scalars at x' . 

To carry out point separation on (5.9), we first define the differential operator 

%b = ^ (1 - 20 {ga^'Va'Vh + gb^'VaVb') + (2^ - gabg^'VcVd' 

(VaVb + ga'^'ob'^'^a'^b') + Cgab (V^V^ + Vc'V"') 

H [Rab - ^gabRj - ^m^gab (5.10) 
from which we obtain the classical stress tensor as 

Tabix) = lim r„b0(x)0(x'). (5.11) 

That the classical tensor field no longer appears as a product of scalar fields at a single point allows 

a smooth transition to the quantum tensor field. From the viewpoint of the stress tensor, the 
separation of points is an artificial construct so when promoting the classical field to a quantum one, 
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neither point should be favored. The product of field configurations is taken to be the symmetrized 
operator product, denoted by curly brackets: 

<t>{x)4>{y) ^ \ {m,4>{y)] = \ {kx)ky) + ky)kx)) (5.i2) 

With this, the point separated stress energy tensor operator is defined as 

f„6(x,x') = \%b[^{x)A{x')] . (5.13) 

While the classical stress tensor was defined at the coincidence limit x' x, we cannot attach any 
physical meaning to the quantum stress tensor at one point until the issue of regularization is dealt 
with, which will happen in the next section. For now, we will maintain point separation so as to 
have a mathematically meaningful operator. 

The expectation value of the point-separated stress tensor can now be taken. This amounts 
to replacing the field operators by their expectation value, which is given by the Hadamard (or 
Schwinger) function 

GW(x,x') = ({to,<^(x')}). (5.14) 
and the point-separated stress tensor is defined as 

{fai,{x,x')) = ^TabG^^\x,x') (5.15) 

where, since 7^^ is a differential operator, it can be taken "outside" the expectation value. The 
expectation value of the point-separated quantum stress tensor for a free, massless (m = 0) con- 
formally coupled (^ = 1/6) scalar field on a four dimension spacetime with scalar curvature R 
is 

(^b-^^ffaft) (5.16) 

5.2.1 Finiteness of Noise Kernel 

We now turn our attention to the noise kernel introduced in Eq. (3.11), which is the symmetrized 
product of the (mean subtracted) stress tensor operator: 

8Nab,c'd'{x,y) = {[fab{x) - {fab{x)),fc'd'{y) - {fc'd'iy))]) 

= {{fM,t,d'{y)})-'2{fab{x)){t,d'{y)) (5.17) 

Since Tab{x) defined at one point can be ill-behaved as it is generally divergent, one can question 
the soundness of these quantities. But as will be shown later, the noise kernel is finite for y ^ x. All 
field operator products present in the first expectation value that could be divergent are canceled 
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by similar products in the second term. We will replace each of the stress tensor operators in the 
above expression for the noise kernel by their point separated versions, effectively separating the two 
points {x, y) into the four points {x, x', y, y'). This will allow us to express the noise kernel in terms 
of a pair of differential operators acting on a combination of four and two point functions. Wick's 
theorem will allow the four point functions to be re-expressed in terms of two point functions. Prom 
this we see that all possible divergences for y ^ x will cancel. When the coincidence limit is taken 
divergences do occur. The above procedure will allow us to isolate the divergences and obtain a 
finite result. 

Taking the point-separated quantities as more basic, one should replace each of the stress tensor 
operators in the above with the corresponding point separated version (5.13), with 7^^ acting at x 
and x' and T^id' acting at y and y' . In this framework the noise kernel is defined as 

8Nab,c'd'{x,y) = lim lim Tab%'d' G{x,x',y,y') (5.18) 

where the four point function is 

G{x,x',y,y') = ^[{{{4>{x)J{x')} , {Hy),ky')}}) 

-2{{kx),kx')}){{ky),ky')})] . (5.19) 

We assume the pairs {x,x') and {y,y') are each within their respective Riemann normal coordinate 
neighborhoods so as to avoid problems that possible geodesic caustics might be present. When wc 
later turn our attention to computing the limit y ^ x, after issues of regularization arc addressed, we 
will want to assume all four points are within the same Riemann normal coordinate neighborhood. 

Wick's theorem, for the case of free fields which we are considering, gives the simple product 
four point function in terms of a sum of products of Wightman functions (we use the shorthand 
notation G^y = G+{x,y) = {4>{x) 4>{y)))- 

{4>{x) 4>{y) 4>{x') ${y')) = Gxy' Gy^' + G^x' Gyy' + G^y G^'y' (5.20) 

Expanding out the anti-commutators in (5.19) and applying Wick's theorem, the four point function 
becomes 

G{x, x' , y, y') = Gxy' Gx'y + Gxy Gx'y' + Gyx' Gy'x + Gyx Gy'x'- (5-21) 

We can now easily see that the noise kernel defined via this function is indeed well defined for the 
limit {x',y') {x,y): 

G{x,x,y,y) = 2(Gly + Glx). (5.22) 

From this we can see that the noise kernel is also well defined for y ^ x; any divergence present 
in the first expectation value of (5.19) have been cancelled by those present in the pair of Green 
functions in the second term, in agreement with the results of section 3. 

5.2.2 Explicit Form of the Noise Kernel 

We will let the points separated for a while so we can keep track of which covariant derivative 
acts on which arguments of which Wightman function. As an example (the complete calculation is 
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quite long), consider the result of the first set of covariant derivative operators in the differential 
operator (5.10), from both Tab and Tc'd', acting on G{x,x' ,y,y'): 

i (1 - 2Cf [ga^'Vp^Vb + 56^"Vp.Va) 

X (5c'^"'Vg- Vd/ + <7d'^"'Vg- Vc') G{x, x', y, y'). (5.23) 

(Our notation is that Va acts at x, Vc' at y, Vf," at x' and V d'" at y'). Expanding out the 
differential operator above, we can determine which derivatives act on which Wightman function: 



X 



5c'^ ff^ a {Gxy';bp"' Gx'y;q"d' + Gxy;bd' G; 



'^i' ■ n" rt'" 



x'y' ■,q"p 



+ Gyx' 


■,q"d' Gyix;bp" 


+ Gyx;bd' Gy'x' 


■qiipiii) 


a {Gxy' -.bp'" Gx'y;q"d 


+ Gxy^bc' Gx'y' 


q"p"' 


+ Gyx' 


■,q"c' Gyix-bp'" 


"1" GyXjbc' Gylx' 


■,q"p"') 


b {Gxy' 


ap'" Gx'y;q"d' 


~l~ Gxy\ad' Gx'y' 


■,q"p"' 


+ Gyx' 


\q" d' Gyix;ap" 


' ~l~ Gyx-^ad' Gylx 


'■,q"p"' 


b {Gxy' 


ap'" Gx'y;q" c! 


~l" Gxy-ac' Gx'y' 


■,q"p"' 


+ Gyx' 


•,q" c'Gy'x;ap"' 


'^Gyx-ac'Gy'x';q 


"p"')] 



(5.24) 

If we now let x' ^ x and y' ^ y the contribution to the noise kernel is (including the factor of | 
present in the definition of the noise kernel): 

g |(1 ~ {Gxy-ad' Gxy;bc' + Gxy;ac' Gxy-bd') 

+ (1 ~ {Gyx;ad' Gyx;bc' + Gyx-ac' Gyx;bd')^ ■ (5.25) 

That this term can be written as the sum of a part involving Gxy and one involving Gyx is a general 
property of the entire noise kernel. It thus takes the form 

Nabc'd'{x, y) = Nabdd' [G+{x, y)] + Nabdd' [G+{y, x)] . (5.26) 

We will present the form of the functional Nabc'd' [G\ shortly. First we note, for x and y time-like 
separated, the above split of the noise kernel allows us to express it in terms of the Feynman (time 
ordered) Green function GF{x,y) and the Dyson (anti-time ordered) Green function GD{x,y). 

Nabc'd' {x, y) = Nabc'd' [Gf{x, y)] + Nabc'd' [Gd{x, y)] . (5.27) 

This can be connected with the zeta function approach to this problem [166] as follows: Recall 
when the quantum stress tensor fluctuations determined in the Euclidean section is analytically 
continued back to Lorentzian signature (r it), the time ordered product results. On the other 
hand, if the continuation is r — > —it, the anti-time ordered product results. With this in mind, the 
noise kernel is seen to be related to the quantum stress tensor fluctuations derived via the effective 
action as 



IGNabc'd' = AT^bc'd' . . . + ^T^bc'd' (5-28) 



2 



t=—iT,t'=—iT' 

The complete form of the functional Nabc'd' [G] is 



t=iT,t'=iT' 



Nabc'd' [G] = Nabc'd' [G] + gabNc'd' [G] + Qc'd'Kb [G] + QabQc'd'N [G] . (5.29) 
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with 



8Nabc'd' [G] = {^ — '^0'^{G.c'bG;d'a + G-c'aG-d'b) 
+4^^ {G-c'd' G-ab + GG-ab&d') 

—2^{1 — 2^){G-bG-c' ad' +G-aG-c'bd'+G-d'G-abc'+G-c'G-abd') 

+2e(l-20(G;a G.bRc'd' + G.c' G,d' Rab) 

—'^^{G-ab Rc'd' + G-c'd' Rab) G + 2^^ Rc'd' RabG^ , 



(5.30) 



8N',,[G] =2(1-20 



-4e 



2^--'j G-^ G-abp'+^ (G-p'^ G-ab + GG-abp'^ 



+2e 



+ ei?') [(1 - 2 - 2 G e G,ab] 

G.p,Gy +2G^G.y' 



2i- 



Rab 



-{rr? + iR!)iRabG\ 



(5.31) 



SN [G] = 2 (^2^-0 G.,p,, Gf^ + (G;/ G,^^^G G.J/) 
+4^(2^-^) (G,pG,/'^ + GP'G,,\) 

? + ^i?) + (m^ + ei?') G.p G'P 



-(2^-2 



-2e [(m' + Cr) Gy' + (m^ + ^i?') G;/ 
+i(m2+^i?)(m2 + Ci?')G^. 



(5.32) 



5.2.3 Trace of the Noise Kernel 



One of the most interesting and surprising results to come out of the investigations undertaken in 
the 1970's of the quantum stress tensor was the discovery of the trace anomaly [217, 218]. When 
the trace of the stress tensor T = g^'^T^b is evaluated for a field configuration that satisties the 
field equation (3.2) the trace is seen to vanish for massless conformally coupled fields. When this 
analysis is carried over to the renormalized expectation value of the quantum stress tensor, the trace 
no longer vanishes. Wald [153] showed this was due to the failure of the renormalized Hadamard 
function Gren(a^, a^') to be symmetric in x and x', implying it does not necessarily satisfy the field 
equation (3.2) in the variable x' . The definition of Gj:en{x,x') in the context of point separation 
will come next.) 

With this in mind, we can now determine the noise associated with the trace. Taking the trace 
at both points x and y of the noise kernel functional (5.27): 

iV[G]=5"'/'*'A^a6c'd'[G'] 



= -3G^ 

9^' 



+ i^i?) G.y' + (m^ + i^i^') G.J 



G.y G,pP + GG,p%:P 



1 



+ [m^ + -^R\ [m^ + -^R!\G 
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+3Q-e) [3(^l-^^G.ypGyP-3^(G.pGyp' + G,jy G,/^'] 



+ U' + Ur) G.yGy + (m' + UR') G.,G'P 



For the massless conformal case, this reduces to 

[G] = {rR'G^ - 6G (Ra' + R'a) G + 18 [(DG) (d'G) + □'□G] }, 



(5.33) 



(5.34) 



which holds for any function G{x, y). For G being the Green function, it satisfies the field equation 
(3.2): 

UG = {rr? + iR)G. (5.35) 

Wc will only assume the Green function satisfies the field equation in its first variable. Using 
the fact W R = (because the covariant derivatives act at a different point than at which R is 
supported), it follows that 

U'UG = {m^ + iR)u'G. (5.36) 
With these results, the noise kernel trace becomes 



N\G\ 



(i-3e) + 3ii (^-e) e 



X [g" (2m2 + i?'c) + (l-60 G,^,G.f -^GiG,/ 
+^ - e) [3 (2 + R! ^ G,p G'P - 18 ^ G.p G.y^' 



+18 



(5.37) 



which vanishes for the massless conformal case. We have thus shown, based solely on the definition 
of the point separated noise kernel, there is no noise associated with the trace anomaly. This 
result obtained in Ref. [8] is completely general since it is assumed that the Green function is only 
satisfying the field equations in its first variable; an alternative proof of this result was given in Ref. 
[6]. This condition holds not just for the classical field case, but also for the regularized quantum 
case, where one does not expect the Green function to satisfy the field equation in both variables. 
One can sec this result from the simple observation used in section 3: since the trace anomaly 
is known to be locally determined and quantum state independent, whereas the noise present in 
the quantum field is non-local, it is hard to find a noise associated with it. This general result 
is in agreement with previous findings [87, 85, 89], derived from the Feynman- Vernon influence 
functional formalism [12, 13] for some particular cases. 



6 Metric fluctuations in Minkowski spacetime 

Although the Minkowski vacuum is an eigenstatc of the total four-momentum operator of a field in 
Minkowski spacetime, it is not an eigenstate of the stress-energy operator. Hence, even for those 
solutions of semiclassical gravity such as the Minkowski metric, for which the expectation value of 
the stress-energy operator can always be chosen to be zero, the fluctuations of this operator are 
non-vanishing. This fact leads to consider the stochastic metric perturbations induced by these 
fluctuations. 
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Here we derive the Einstein-Langevin equation for the metric perturbations in a Minkowski 
background. We solve this equation for the hnearized Einstein tensor and compute the associated 
two-point correlation functions, as well as, the two-point correlation functions for the metric pertur- 
bations. Even though, in this case, we expect to have negligibly small values for these correlation 
functions for points separated by lengths larger than the Planck length, there are several reasons 
why it is worth carrying out this calculation. 

On the one hand, these are the first back-reaction solutions of the full Einstein-Langevin equa- 
tion. There are analogous solutions to a "reduced" version of this equation inspired in a "mini- 
superspace" model [90, 91], and there is also a previous attempt to obtain a solution to the Einstein- 
Langevin equation in Ref. [89], but, there, the non-local terms in the Einstein-Langevin equation 
were neglected. 

On the other hand, the results of this calculation, which confirm our expectations that gravita- 
tional fluctuations are negligible at length scales larger than the Planck length, but also predict that 
the fluctuations are strongly suppressed on small scales, can be considered a first test of stochastic 
semiclassical gravity. In addition, these results reveal an important connection between stochastic 
gravity and the large N expansion of quantum gravity. We can also extract conclusions on the pos- 
sible qualitative behavior of the solutions to the Einstein-Langevin equation. Thus, it is interesting 
to note that the correlation functions at short scales are characterized by correlation lengths of the 
order of the Planck length; furthermore, such correlation lengths enter in a non-analytic way in the 
correlation functions. 

We advise the reader that his section is rather technical since it deals with an explicit non trivial 
back reaction computation in stochastic gravity. We tried to make it reasonable self-contained and 
detailed, however a more detailed exposition can be found in Ref. [94]. 

6.1 Perturbations around Minkowski spacetime 

The Minkowski metric rjab, in a manifold M. which is topologically M^, and the usual Minkowski 
vacuum, denoted as jO), are the class of simplest solutions to the semiclassical Einstein equation 
(3.7), the so called trivial solutions of semiclassical gravity [219]. They constitute the ground state 
of semiclassical gravity. In fact, we can always choose a renormalization scheme in which the 
renormalized expectation value (0| Tg'' [ry]|0) = 0. Thus, Minkowski spacetime (IR'',77a(,) and the 
vacuum state |0) are a solution to the semiclassical Einstein equation with renormalized cosmological 
constant A = 0. The fact that the vacuum expectation value of the renormalized stress-energy 
operator in Minkowski spacetime should vanish was originally proposed by Wald [152] and it may 
be understood as a renormalization convention [33, 35]. Note that other possible solutions of 
semiclassical gravity with zero vacuum expectation value of the stress-energy tensor are the exact 
gravitational plane waves, since they are known to be vacuum solutions of Einstein equations which 
induce neither particle creation nor vacuum polarization [220, 221, 222]. 

As we have already mentioned the vacuum |0) is an eigenstate of the total four- momentum 
operator in Minkowski spacetime, but not an eigenstate of T^[ri]. Hence, even in the Minkowski 
background, there are quantum fluctuations in the stress-energy tensor and, as a result, the noise 
kernel does not vanish. This fact leads to consider the stochastic corrections to this class of trivial 
solutions of semiclassical gravity. Since, in this case, the Wightman and Feynman functions (4.23), 
their values in the two-point coincidence limit, and the products of derivatives of two of such 
functions appearing in expressions (4.24) and (4.25) are known in dimensional regularization, we 
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can compute the Einstein-Langevin equation using the methods outhned in sections 3 and 4. 

To perform explicit calculations it is convenient to work in a global inertial coordinate system 
{x^^} and in the associated basis, in which the components of the flat metric are simply ry^^ = 
diag(— 1, 1, . . . , 1). In Minkowski spacetime, the components of the classical stress-energy tensor 
(3.3) reduce to 

T^^'lri, 0] = - ^ r]'"'df^dp(l) - ^ 77'^^ + C iv^''^ - d^d") 4>^ , (6.1) 

where = 8^0'^, and the formal expression for the components of the corresponding "operator" in 
dimensional regularization, see Eq. (3.4), is 

rr W = lid'k, d^k) + T^'^^l (6.2) 

where Vf^^ is the differential operator (3.5), with g^i, = rj/j,^, R^i, = 0, and V/j, = d^. The field 

4>nix) is the field operator in the Heisenberg representation in a n-dimensional Minkowski spacetime, 
which satisfies the Klein-Gordon equation (3.2). We use here a stress-energy tensor which differs 
from the canonical one, which corresponds to ^ = 0, both tensors, however, define the same total 
momentum. 

The Wightman and Feynman functions (4.23) when g^i, = 7)^1,, are well known: 

G+(x,y) = iA+ix - y), G^^{x,y) = A^Jx - y), (6.3) 

with 



/a k 



(2vr) 

where /c^ = "q^j^yk^k^ and kx = rjfj^yk'^x'^ . Note that the derivatives of these functions satisfy 
5^A+(x — y) = 9^A+(a; — y) and 9-^A+(x — y) = — 5^A+(x — y), and similarly for the Feynman 
propagator ApJ^x — y). 

To write down the semiclassical Einstein equation (3.7) in n-dimensions for this case, we need 
to compute the vacuum expectation value of the stress-energy operator components (6.2). Since, 
from (6.3), we have that (0|^^(x)|0) = iAj,^(0) = iA+(0), which is a constant (independent of x), 
we have simply 

io\Trim=-^Jj^ ,,_^„,,_,^ = '-[-^) r(--), (6.5) 

where the integrals in dimensional regularization have been computed in the standard way (see 
Ref. [94]) and where r(z) is the Euler's gamma function. The semiclassical Einstein equation (3.7) 
in n-dimensions before renormalization reduces now to 

A« r^^^^ = f,-in-^){0\fli^[rj]\o). (6.6) 
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This equation, thus, simply sets the value of the bare coupling constant Kb/Gb- Note, from 
(6.5), that in order to have (0| Tj^^ |0)[r/]=0. the rcnornializcd and regularized stress-energy tensor 
"operator" for a scalar field in Minkowski spacetime, sec Eq. (3.6), has to be defined as 

which corresponds to a renormalization of the cosmological constant 

As A 2 



Gb G TT n(n-2) 

where 



Kn + 0{n-4), (6.8) 



being 7 the Euler's constant. In the case of a massless scalar field, = 0, one simply has 
^b/Gb = A/G. Introducing this renormalized coupling constant into Eq. (6.6), we can take the 
limit n— >-4. We find that, for (IR^,r/a6, |0)) to satisfy the semiclassical Einstein equation, we must 
take A = 0. 

We can now write down the Einstein-Langevin equations for the components h^i, of the stochas- 
tic metric perturbation in dimensional regularization. In our case, using (0|(j!)^(a;)|0) = iAji^^(O) and 
the explicit expression of Eq. (4.20) we obtain 



1 



- asL>(i^'^'^(x) - 2^bB(i)'*'^(x) 
3 



SttGs 

-^&'^i'-{x)ii-^^-^\/:^^^{<d) + ]^jd^yix-^^-^^^^^ (6.10) 

The indices in h^y are raised with the Minkowski metric and h = h^,, and here a superindex (1) 
denotes the components of a tensor linearized around the fiat metric. Note that in n-dimensions 
the two-point correlation functions for the field ^'^'^ is written as 

{e''{x)e^{y))s = M-'^"-'^iVr"''(x,y), (6.11) 
Explicit expressions for and i?'^'^*^ are given by 

D^'^'^ix) = \j^r^Kp{x\ B^'^'^'^ix) = 2J^^^J^fKp{x), (6.12) 

with the differential operators T^^" = rj'^'n^ - d^d"^ and J^"'^ = 3J"^^"jf - T!^"" . 



6.2 The kernels in the Minkowski background 

Since the two kernels (4.22) are free of ultraviolet divergences in the limit n — > 4, we can deal 
directly with the F^'^'^l^ix - y) = lim„^4 /x'^ ("-^) F'^""''^ in Eq. (4.21). The kernels N'^'"^^{x,y) = 
Ke F^^*^^ (x — y) and H^^'^^{x,y) = Im F^'^°'l^ (x — y) are actually the components of the "physical" 
noise and dissipation kernels that will appear in the Einstein-Langevin equations once the renor- 
malization procedure has been carried out. The bi-tensor F^^^°^l^ can be expressed in terms of the 
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Wightman function in four spacetime dimensions, according to (4.24). The different terms in this 
kernel can be easily computed using the integrals 



lip) = I j—^ 5{k^ + m^) e{-k^) S[ik - pf + 0(A;O - (6.13) 



and //^i-/^''(p) which are defined as the previous one by inserting the momenta A;^^ . . .k^'^ with 
r=l,2,3,4 in the integrand. All these integral can be expressed in terms of /(p); see Ref. [94] for 
the explicit expressions. It is convenient to separate /(p) in its even and odd parts with respect to 
the variables p^ as 

^(p)=4(p)+/a(p), (6.14) 
where Is{—p) = Is{p) and Ip,{—p) = —Ia{p)- These two functions are explicitly given by 



= ^( V - 4m2) J 1 + 4 



8(27r)3 ^ ^ ^ V p 



2 



Ia{p) = g^2^ sign/ e{-p^ - An?) ^ 1 + 4 (6.15) 



After some manipulations, we find 



F^^'^^i.) = g J 7^^ e-- (l + 4 ^) ' Hp) 



(27r)4 \^ p^ 

+ ^ ^re'^ / e-- (3 A^ + ^J lip), (6.16) 

where = ^ — 1/6. The real and imaginary parts of the last expression, which yield the noise 
and dissipation kernels, are easily recognized as the terms containing Isip) and Ia{p), respectively. 
To write them explicitly, it is useful to introduce the new kernels 



m? ( . ^ m? 



2 



-i f (fp 0/ J2 A^2\ 



DAix; m^)^^J-^ e^f- signp" ^( V - 4.m' 

2 



D^ix- m\AO = T^J-0^ e'P- signp^ 0(-/ _ 4m^) 



(2 

x/l + 4^(3AC + ^y (6-17) 
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and we finally get 

N^^^^{x,y) = ^JT"''iVA(x-y;m2) + J•^^,"'3iVB(a;-y;m^ AO, 

Hr^{x, y) = ^J^r^DA{x-y; m^) + J^^'^J^fD^ix-y; m\ A^- (6.18) 

Notice that the noise and dissipation kernels defined in (6.17) are actually real because, for the 
noise kernels, only the cos pa; terms of the exponentials e*^^ contribute to the integrals, and, for the 
dissipation kernels, the only contribution of such exponentials comes from the i sin pa; terms. 

The evaluation of the kernel H^'^"'^ {x,y) is a more involved task. Since this kernel contains 
divergences in the limit n— >4, we use dimensional regularization. Using Eq. (4.25), this kernel can 
be written in terms of the Feynman propagator (6.4) as 



where 



^_(„_4)^^/.^.a/3^^^ y) = Imi^/^'^"^(x - y), (6.19) 



X'*'^°^(a;) = -//-("-^)|2a'^a("A^„(a;)5^)a'^A^„(a;) 
+2P^''(a"A^„(x)a'5A^„(a;)) + 2P"^(a^A^„(x) a'^A^„(a 



+2Vi"'V'^^ (a|^(x)) + rii"'d^'^Aj.^{x)d^^ +ri'^^d^''A^^{x)d''^ 

+A,„(0)(,?^'^P"^+<'^P^'^) + ir?'^V''(A.„(a;)D-m2A,JO))]5-(x)|. (6.20) 
Let us define the integrals 

J„(p) = 11 ^^jj^ (P + m2-ze) p-p)2 + m2-k]' ^^'^^^ 

and J^i -'^''(p) obtained by inserting the momenta fc'^i . . . k^'' into the previous integral, together 
with 

(If (FT^' 

and Iq^'"^"" which are also obtained by inserting momenta in the integrand. Then, the different 
terms in Eq. (6.20) can be computed. These integrals are explicitly given in Ref. [94]. It is found 
that Iq^^ = and the remaining integrals can be written in terms of /o„ and Jn{p)- It is useful to 
introduce the projector P^^ orthogonal to p^ and the tensor pt^^'^l^ as 

p2 pnv ^ _ plJ'p'^ ^ pnvajd ^ 3p/u(a pjd^ _ p^iv paji ^ (6.23) 

then the action of the operator J^i^^ is simply written as J^^'^ J d^pe'^^ /(p) = — / d^pe*^^ f{p)p^P'^'^ 
where /(p) is an arbitrary function of p'^. 

Finally after a rather long but straightforward calculation, and after expanding around n = 4, 
we get, 



(47r)2 



90 
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+ 



2 m2 



3 (n-2) 



n(n — 2) 



180^ 7 (27r)'^ 



e*p- i + 4::_ 



-m 



2 1 
135 ^ 27 ^ ^ 



A„(x)'.+0(n-4), 



where has been defined in (6.9), and and A„(x) are given by 

n2 



(6.24) 



(6.25) 



(6.26) 



where ip{p'^) = Jq^ da In |1 + ^ a(l— Q;)|. The imaginary part of (6.24) gives the kernel components 

pi^-i^-^) H^^"^!^ {x,y), according to (6.19). It can be easily obtained multiplying this expression by 
—i and retaining only the real part, <^(p^), of the function 0(p^). 

6.3 The Einstein-Langevin equation 

With the previous results for the kernels we can now write the n-dimensional Einstein-Langevin 
equation (6.10), previous to the renormalization. Taking also into account Eqs. (6.5) and (6.6), we 
may finally write: 

^_G(^)^'^(a;) - ^ aBD^''^''{x) - 2/?Bi?(^)^-(x) 



+ 



+ 



(4^)2 

1 



^ (n-2) 90 ^ 



(x) 



16 



28807r2 15 



(Pp 
(27r)" 



. m 



1 + 4— D^'^^^^iy) 



p^ 



+ 10(^3A^+-^J B«'^'^(y) 
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(6.27) 



Notice that the terms containing tlie bare cosmological constant have canceUed. These equations 
can now be renormaUzed, that is, we can now write the bare couphng constants as renormahzed 
couphng constants plus some suitably chosen counterterms and take the limit n — ^ 4. In order to 
carry out such a procedure, it is convenient to distinguish between massive and massless scalar 
fields. The details of the calculation can be found in Ref. [94]. 
It is convenient to introduce the two new kernels 



HA{x;m^) 



1 



4807r2y (27r)4 
ZTTsignp" 9 ( —p^ — 4m^ ) 



1 + 4- 



m 



p^ 



'1 + 4- 



8 

3^ 



HB{x;m'',AO = 



d^p 

727r2j (27r)4^ 



3AC + - 



m 



— iTTsign p° 9 ( -p"^ -4m'^)\l l+4^+ip{p 



p^ 



Im 
6 p2 



(6.28) 



where ^pip^) is given by the restriction to n = 4 of expression (6.25). The renormalized coupling 
constants 1/G, a and (5 are easily computed as it was done in Eq. (6.8). Substituting their 
expressions into Eq. (6.27), we can take the limit n— >4, using the fact that, for n = 4, W^^^^^^x) = 
(3/2) ^'^'^'^(x), we obtain the corresponding semiclassical Einstein-Langevin equation. 

For the massless case one needs the limit m^O of equation (6.27). In this case it is convenient 
to separate Kn in (6.9) as = + ^ \n{m? / jx^) + 0(n— 4), where 



Kn. — 



ei \ n-4 



(n— 4) v47r 
and use that, from Eq. (6.25), we have 



lim 



(/?(p^) + In(mVAt^) 



V47r 



-2 + In 



+ 0(n-4), 



(6.29) 



(6.30) 



The coupling constants are then easily renormalized. We note that in the massless limit, the 
Newtonian gravitational constant is not renormalized and, in the conformal coupling case, = 0, 
we have that Pb = P- Note also that, by making m = in (6.17), the noise and dissipation kernels 
can be written as 



DA{x;m^ = 0) = D{x), 



NB{x;m'^ = 0,AC) = 60AfN{x), 
L>B(a;;m2 = 0,AO = 60AfD{x), 



where 



(6.31) 
(6.32) 
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It is also convenient to introduce the new kernel 

1 f d^p 



This kernel is real and can be written as the sum of an even part and an odd part in the variables 
x'^, where the odd part is the dissipation kernel D{x). The Fourier transforms (6.32) and (6.33) 
can actually be computed and, thus, in this case we have explicit expressions for the kernels in 
position space; see, for instance, Refs. [223, 224, 225]. 

Finally, the Einstein-Langevin equation for the physical stochastic perturbations h^y can be 
written in both cases, for m^^O and for m = 0, as 

^ G'^'^'^(x) - 2 {aA^''^i'''{x) + ^5(1)'^^ (x)) 



SttG 



+\ j [H^{x-y)A^'^i'^{y) + H^{x-y)B('^>'%y)\ = e^ix), (6.34) 

where in terms of the renormalized constants a and (3 the new constants are a = a + (36007r^)~^ 
and (3 = P - (1/12 - 5AO(28807r2)-i. The kernels Ha{x) and Hb{x) are given by Eqs. (6.28) 
when m 7^ 0, and Ha{x) = H{x;^'^), Hb{x) = 60 A^^i?(a;; //^) when m = 0. In the massless 
case, we can use the arbitrariness of the mass scale fi to eliminate one of the parameters a or f3. 
The components of the Gaussian stochastic source ^^'^ have zero mean value and their two-point 
correlation functions are given by {^'"'{x)^'^f^{y))s = N'""^f^{x,y), where the noise kernel is given 
in (6.18), which in the massless case reduces to (6.31). 

It is interesting to consider the massless conformally coupled scalar field, i.e., the case A^ = 0, of 
particular interest because of its similarities with the electromagnetic field, and also because of its 
interest in cosmology: massive fields become conformally invariant when their masses arc negligible 
compared to the spacetime curvature. We have already mentioned that for a conformally coupled, 
field, the stochastic source tensor must be traceless (up to first order in perturbation theory around 
semiclassical gravity), in the sense that the stochastic variable i'^ = 'n^u^^'^ behaves deterministically 
as a vanishing scalar field. This can be directly checked by noticing, from Eqs. (6.18) and (6.31), 
that, when A^ = 0, one has {i'^{x)e,'^'^ {y)) s = 0, since J^/^ = 3D and J^i^'^T^ = UJ^'^P . The Einstein- 
Langevin equations for this particular case (and generalized to a spatially flat Robertson- Walker 
background) were first obtained in Ref. [89], where the coupling constant (3 was fixed to be zero. See 
also Ref. [4] for a discussion of this result and its connection to the problem of structure formation 
in the trace anomaly driven inflation [113, 114, 115]. 

Note that the expectation value of the renormalized stress-energy tensor for a scalar field can be 
obtained by comparing Eq. (6.34) with the Einstein-Langevin equation (3.14), its explicit expression 
is given in Ref. [94]. The results agree with the general form found by Horowitz [225, 226] using an 
axiomatic approach and coincides with that given in Ref. [219]. The particular cases of conformal 
coupling, A^ = 0, and minimal coupling, A^ = — 1/6, are also in agreement with the results for this 
cases given in Refs. [225, 226, 227, 84, 24], modulo local terms proportional to A^'^^^^^ and B^^^^^^ 
due to different choices of the renormalization scheme. For the case of a massive minimally coupled 
scalar field, A^ = — 1/6, our result is equivalent to that of Ref. [228]. 
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6.4 Correlation functions for gravitational perturbations 

Here we solve the Einstein-Langevin equations (6.34) for the components G^^"*^" of the linearized 
Einstein tensor. Then we use these solutions to compute the corresponding two-point correlation 
functions, which give a measure of the gravitational fluctuations predicted by the stochastic semi- 
classical theory of gravity in the present case. Since the linearized Einstein tensor is invariant 
under gauge transformations of the metric perturbations, these two-point correlation functions are 
also gauge invariant. Once we have computed the two-point correlation functions for the linearized 
Einstein tensor, we find the solutions for the metric perturbations and compute the associated 
two-point correlation functions. The procedure used to solve the Einstein-Langevin equation is 
similar to the one used by Horowitz [225], see also Ref. [219], to analyze the stability of Minkowski 
spacctimc in scmiclassical gravity. 

We first note that the tensors A'^^^^^ and B^^^^^^ can be written in terms of G'^^^^^ as 

^(DM^^ = I (jrM'^G^i)^ - J^^G^'^^""), B^'^t"" = 2J^^'^G<^'2, (6.35) 

where we have used that 3n = J^^. Therefore, the Einstein-Langevin equation (6.34) can be seen as 
a linear intcgro-differential stochastic equation for the components G*^'^*^. In order to find solutions 
to Eq. (6.34), it is convenient to Fourier transform. With the convention f{p) = J d'^xe~^^^ f{x) 
for a given field f{x), one finds, from (6.35), 

A^'^^"' {p) = 2p'^&''>^"' (p) - ^^P^''&''> I {p) , 

B^^)!''' (p) = -2/p/^'^G(i) ^ {p) . (6.36) 
The Fourier transform of the Einstein-Langevin Eq. (6.34) now reads 

^'^ a/3(p) G^'^'^ip) = STrGTip), (6.37) 

where 

F'^apiP) ^ Mp) S^J'p) + F2{p)p'P'^''r,^^, (6.38) 



with 

Fi{p) = 1 + 167rGp2 



^HAip) - 2a 



F2(p) = -y7rG 



(6.39) 



In the Fourier transformed Einstein-Langevin Eq. (6.37), ^'^'^{p), the Fourier transform of ^'^"{x), 
is a Gaussian stochastic source of zero average and 

{e''{p)eHp'))s = (2'K)U^{j>+p')m^-^{j>), (6.40) 

where we have introduced the Fourier transform of the noise kernel. The explicit expression for 
j^nua/B^p^j is found from (6.18) and (6.17) to be 

I / o^2 

N'^'^'^ip) = ^Gi-P' 




(6.41) 
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which in the massless case reduces to 

hm m^'^i^ip) = ^— ei-p"^ 

4807r ^ 



(6.42) 



6.4.1 Correlation functions for the linearized Einstein tensor 

In general, we can write G^^^^''' = {G^^'>i"')s + Gf'""", where Of^'^" is a solution to Eqs. (6.34) 
with zero average, or (6.37) in the Fourier transformed version. The averages {G'-'^^f^^)s must be 
a solution of the linearized semiclassical Einstein equations obtained by averaging Eqs. (6.34), or 
(6.37). Solutions to these equations (specially in the massless case, m = 0) have been studied by 
several authors [225, 226, 229, 230, 231, 232, 233, 181, 234, 24, 219], particularly in connection with 
the problem of the stability of the ground state of semiclassical gravity. The two-point correlation 
functions for the linearized Einstein tensor are defined by 



= (G«'^''(x)G«"^(a;')).. 



(6.43) 



Now we shall seek the family of solutions to the Einstein-Langevin equation which can be 
written as a linear functional of the stochastic source and whose Fourier transform, G'-^^^'^{p), 
depends locally on ^"'^(p). Each of such solutions is a Gaussian stochastic field and, thus, it can be 
completely characterized by the averages (G'^'^'')^ and the two-point correlation functions (6.43). 
For such a family of solutions, Gf-'^'^{p) is the most general solution to Eq. (6.37) which is linear, 
homogeneous and local in ^"^(p). It can be written as 



Gf'^''{p) = 87rGD^\pip)eP{p), 



(6.44) 



apip) components of a Lorentz invariant tensor field distribution in Minkowski 



where D^'^ 

spacetime (by "Lorentz invariant" we mean invariant under the transformations of the orthochronous 
Lorentz subgroup; see Ref. [225] for more details on the definition and properties of these tensor 
distributions), symmetric under the interchanges a ^ (3 and /x <-> z^, which is the most general 
solution of 



(6.45) 



In addition, we must impose the conservation condition: PvGf^^^ {p) = 0, where this zero must be 
understood as a stochastic variable which behaves deterministically as a zero vector field. We can 



write D 



al3 



(P) 



ip) + 



ip) 



(p), where D'^'^ ^p{p) is a particular solution to Eq. (6.45) and 
to th( 

G^p^'^'^ip) + &i^^^'^{p). To find the particular solution, we try an ansatz of 



apip) most general solution to the homogeneous equation. Consequently, see Eq. (6.44), 

we can write G^^"^'^ 
the form 



D^\pip) = diip) 5[J^0) + d2{p)p^P^''vo.p. (6.46) 
Substituting this ansatz into Eqs. (6.45), it is easy to see that it solves these equations if we take 



di{p) 



Flip) 



d2{p) 



F2{p) 



Fi(jp)F^{jp) 



with 



F^(jp) = Fi{p) + 3p'^F2{p) = 1 - 487rGp2 



\hb{p)-2P 



(6.47) 



(6.48) 
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and where the notation [ ]r means that the zeros of the denominators are regulated with appro- 
priate prescriptions in such a way that di{p) and d2{p) are well defined Lorentz invariant scalar 
distributions. This yields a particular solution to the Einstein-Langevin equations: 



G 



(6.49) 



which, since the stochastic source is conserved, satisfies the conservation condition. Note that, in 
the case of a massless scalar field, m = 0, the above solution has a functional form analogous to 
that of the solutions of linearized semiclassical gravity found in the Appendix of Ref. [219]. Notice 
also that, for a massless conformally coupled field, m = and A^ = 0, the second term on the right 
hand side of Eq. (6.46) will not contribute in the correlation functions (6.43), since in this case the 
stochastic source is traccless. 

A detailed analysis given in Ref. [94] concludes that the homogeneous solution &^^^^ {p) gives no 
contribution to the correlation functions (6.43). Consequently g^''"^l^{x, x') = (G«''''(x)G«"^(x'))s, 
where G^^'^'^{x) is the inverse Fourier transform of (6.49), and the correlation functions (6.43) are 

{Gl'^^^'^ip)&;^''%'))s = 6i2nf^^^^^ (6.50) 

It is easy to see from the above analysis that the prescriptions [ ]r in the factors Dp are irrelevant 
in the last expression and, thus, they can be suppressed. Taking into account that Fi{—p) = F*{p), 
with / = 1,2,3, we get from Eqs. (6.46) and (6.47) 



{&;^^^{p)G^;^-^ip')), 

64(2. f G^ ^-^^^ 
\Fi{p)\ 



N^'^'^'^ip) - p'^P^^'N'^'^PJp) 



Fsip) 



^3*(P) 



I^3(P)| 



%{P) 



(6.51) 



This last expression is well defined as a bi-distribution and can be easily evaluated using Eq. (6.41). 
The final explicit result for the Fourier transformed correlation function for the Einstein tensor is 
thus 



(GWM-(p)(5Wa/3(^'))^ 
5\p+p') 



45^ ' 



\Fi{p)Y 



e{-p^-4m^)\ 1 + A 



p^ 



2\2 -plJ.i'ocP 



+10^3Ae+^j (p'^fpt'^p^f^ 



l-3p- 



Fsip) 



(6.52) 



To obtain the correlation functions in coordinate space, Eq. (6.43), we take the inverse Fourier 
transform. The final result is: 

x') = ^ gA{x " + y ^r^f Gb{x - x'), (6.53) 
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(6.54) 



where Fi{p), I = 1,2,3, are given in (6.39) and (6.48). Notice that, for a massless field (m = 0), we 
have 



Fi{p) = l + 4nGp''H{p;fl^), 



16 

i^2(p) = -y7rG 



(i + 180AO^^(p;A')-6T 



Fsd?) = 1 - 487rG/ 15A^2^(p; /x^) - 2T 



(6.55) 



with i-i ^ exp(19207r^Q;) and T = /? — 60 A^^ a, and where H{p\^'^) is the Fourier transform of 
H{x;fi^) given in (6.33). 



6.4.2 Correlation functions for the metric perturbations 

Starting from the solutions found for the linearized Einstein tensor, which are characterized by 
the two-point correlation functions (6.53) [or, in terms of Fourier transforms, (6.52)], we can now 

solve the equations for the metric perturbations. Working in the harmonic gauge, d^h^'^ = (this 
zero must be understood in a statistical sense) where h^u = hni,—{l/2)r]fj,,^h^, the equations for the 
metric perturbations in terms of the Einstein tensor are 

nh^"'{x) = -2G'-''>^'''{x), (6.56) 

or, in terms of Fourier transforms, p^h>^'^{p) = 2&^^^^'^{p). Similarly to the analysis of the equation 
for the Einstein tensor, we can write h^^ = iji^^) s + h'^^ , where hf^^ is a solution to these equations 
with zero average, and the two-point correlation functions are defined by 

H^''''P{x,x') = {h^''{x)h''^{x'))s - {h'"'{x))s{h''^{x'))s 

= {h^^{x)hf{x'))s. (6.57) 

We seek solutions of the Fourier transform of Eq. (6.56) of the form h'^^ {p) = 2D{p)Gf ^^^ [p), 
where D{p) is a Lorentz invariant scalar distribution in Minkowski spacetime, which is the most 
general solution of p^Dij)) = 1. Note that, since the linearized Einstein tensor is conserved, solutions 
of this form automatically satisfy the harmonic gauge condition. As in the previous subsection, 
we can write D{p) = [l/p'^]r + Dh{p), where Dh{p) is the most general solution to the associated 
homogeneous equation and, correspondingly, we have h^'^ip) = hp^ ip) + {p). However, since 
Dfi{p) has support on the set of points for which p^ = 0, it is easy to see from Eq. (6.52) (from the 
factor 9{—p^ — 4m^)) that (h'^'^ {p)G'^^'^^ {p')) g = and, thus, the two-point correlation functions 
(6.57) can be computed from {Ji'^''{p)hf{p'))s = {hf^''{p)hf{p'))s. From Eq. (6.52) and due to the 
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factor 0{—p^ — 4m^), it is also easy to see that the prescription [ is irrelevant in this correlation 
function and we obtain 



{k''{p)nfip'))s = &;"'^p'))s, (6.58) 

where {&p^^^^{p)&p^°'^{p'))s is given by Eq. (6.52). The right hand side of this equation is a well 
defined bi-distribution, at least for m 7^ (the 9 function provides the suitable cutoff). In the 
massless field case, since the noise kernel is obtained as the limit m — > of the noise kernel for a 
massive field, it seems that the natural prescription to avoid the divergences on the lightcone = 
is a Hadamard finite part, see Refs. [235, 236] for its definition. Taking this prescription, we also 
get a well defined bi-distribution for the massless limit of the last expression. 
The final result for the two-point correlation function for the field h^^ is: 

Anv 327r 

Hi'^'^^ix, x') = ^G^ J^r"'' HAix -x') + ^ ^r-^x^ ^B(ic - x'), (6.59) 

where HAip) = [l/{p'^?] Qa{p) and Hb(p) = [l/(p^)^] Gb(p), with GAip) and GBip) given by (6.54). 
The two-point correlation functions for the metric perturbations can be easily obtained using h^i, = 

V-(V2)?7m-/iS- 

6.4.3 Conformally coupled field 

For a conformally coupled field, i.e., when m = and = 0, the previous correlation functions are 
greatly simplified and can be approximated explicitly in terms of analytic functions. The detailed 
results are given in Ref. [94], here we outline the main features. 

When m = and = we have ^B(a;) =0 and ^a(p) = 0{-p^) \ Fi{p) |~^. Thus the two-point 
correlations functions for the Einstein tensor is 

where H{p,ii'^) = (4807r2)-i \n[-{{jp + ief + p'pi)/ fi"^], see Eq. (6.33). 

To estimate this integral let us consider spacelike separated points {x — x'Y = (0,x — x'), and 
define y = x — x'. We may now formally change the momentum variable p^^ by the dimensionless 
vector s^: p^ = s^/|y|, then the previous integral denominator is |1 -|- 167r(Lp/|y|)^s^i/(s)|^, where 
we have introduced the Planck length Lp = \fG. It is clear that we can consider two regimes: 
(a) when Lp <^ |y|, and (b) when |y| ~ Lp. In case (a) the correlation function, for the 0000 
component, say, will be of the order 

In case (b) when the denominator has zeros a detailed calculation carried out in Ref. [94] shows 
that: 

which indicates an exponential decay at distances around the Planck scale. Thus short scale 
fluctuations are strongly suppressed. 
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For the two-point metric correlation the results are similar. In this case we have 



The integrand has the same behavior of the correlation function of Eq. (6.60) thus matter fields 
tends to suppress the short scale metric perturbations. In this case we find, as for the correlation 
of the Einstein tensor, that for case (a) above we have, 

and for case (b) we have 

It is interesting to write expression (6.61) in an alternative way. If we use the dimensionless 
tensor pt^^°'P introduced in Eq. (6.23), which accounts for the effect of the operator , we can 

write 

H^'^-Hx x') = ^G' M^-^') p,^^P e{-p') 

^ ' ^ 45 J (27r)4 |l + 47rGp2#(p;^2)|2 • ^ " ^ 

This expression allows a direct comparison with the graviton propagator for linearized quantum 
gravity in the expansion found by Tomboulis [182]. One can see that the imaginary part of 
the graviton propagator leads, in fact, to Eq. (6.62). In Ref. [169] it is shown that, in fact, the 
two-point correlation functions for the metric perturbations derived from the Einstein-Langevin 
equation are equivalent to the symmetrized quantum two-point correlation functions for the metric 
fiuctuations in the large N expansion of quantum gravity interacting with N matter fields. 



6.5 Discussion 

The main results of this section are the correlation functions (6.53) and (6.59). In the case of 
a conformal field, the correlation functions of the linearized Einstein tensor have been explicitly 
estimated. From the exponential factors e"'^'/^^ in these results for scales near the Planck length, 
we sec that the correlation functions of the linearized Einstein tensor have the Planck length as 
the correlation length. A similar behavior is found for the correlation functions of the metric 
perturbations. Since these fluctuations are induced by the matter fluctuations we infer that the 
effect of the matter fields is to suppress the fluctuations of the metric at very small scales. On the 
other hand, at scales much larger than the Planck length the induced metric fluctuations are small 
compared with the free graviton propagator which goes like Lp/|y|2, since the action for the free 
graviton goes like Sh ^ J d'^x L'p^hUh 

For background solutions of semiclassical gravity with other scales present apart from the Planck 
scales (for instance, for matter fields in a thermal state), stress-energy fluctuations may be impor- 
tant at larger scales. For such backgrounds, stochastic semiclassical gravity might predict correla- 
tion functions with characteristic correlation lengths larger than the Planck scales. It seems quite 
plausible, nevertheless, that these correlation functions would remain non-analytic in their charac- 
teristic correlation lengths. This would imply that these correlation functions could not be obtained 
from a calculation involving a perturbative expansion in the characteristic correlation lengths. In 
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particular, if these correlation lengths are proportional to the Planck constant h, the gravitational 
correlation functions could not be obtained from an expansion in h. Hence, stochastic semiclassical 
gravity might predict a behavior for gravitational correlation functions different from that of the 
analogous functions in perturbative quantum gravity [198, 199, 200, 201]. This is not necessarily 
inconsistent with having neglected action terms of higher order in ti when considering semiclassical 
gravity as an effective theory [219]. It is, in fact, consistent with the closed connection of stochastic 
gravity with the large N expansion of quantum gravity interacting with N matter fields. 

7 Structure formation 

Cosmological structure formation is a key problem in modern cosmology [237, 238] and inflation 
offers a natural solution to this problem. If an inflationary period is present, the initial seeds for 
the generation of the primordial inhomogeneities that lead to the large scale structure have their 
source in the quantum fluctuations of the inflaton field, the field which is generally responsible for 
driving inflation. Stochastic gravity provides a sound and natural formalism for the derivation of 
the cosmological perturbations generated during inflation. 

In Ref. [112] it was shown that the correlation functions that follow from the Einstein-Langevin 
equation which emerges in the framework of stochastic gravity coincide with that obtained with 
the usual quantization procedures [239] when both the metric perturbations and the inflaton fluc- 
tuations are both linearized. Stochastic gravity, however, can naturally deal with the fluctuations 
of the inflaton field even beyond the linear approximation. 

Here we will illustrate the equivalence with the usual formalism, based on the quantization of 
the linear cosmological and inflaton perturbations, with one of the simplest chaotic inflationary 
models in which the background spacetime is a quasi de Sitter universe [240, 112]. 

7.1 The model 

In this chaotic inflationary model [241] the inflaton field (f) of mass m is described by the following 
Lagrangian density, 

m = ^5"'V„0V6</. + ^m2</,2. (7.1) 

The conditions for the existence of an inflationary period, which is characterized by an accelerated 
cosmological expansion, is that the value of the field over a region with the typical size of the Hubble 
radius is higher than the Planck mass mp. This is because in order to solve the cosmological horizon 
and flatness problem more than 60 e-folds of expansion are needed, to achieve this the scalar field 
should begin with a value higher than 3mp. The inflaton mass is small: as we will see, the large 
scale anisotropics measured in the cosmic background radiation [242] restrict the inflaton mass to 
be of the order of lO^^mp. We will not discuss the naturalness of this inflationary model and we 
will simply assume that if one such region is found (inside a much larger universe) it will inflate to 
become our observable universe. 

Wc want to study the metric perturbations produced by the stress-energy tensor fluctuations 
of the inflaton field on the homogeneous background of a flat Friedmann-Robertson- Walker model, 
described by the cosmological scale factor 0(77), where r] is the conformal time, which is driven by 
the homogeneous infiaton field (f){r]) = {(p). Thus we write the infiaton field in the following form 

= 0(77) + (^(x), (7.2) 
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where 0{x) corresponds to a free massive quantum scalar field with zero expectation value on the 
homogeneous background metric: (i^) = 0. Wc will restrict ourselves to scalar-type metric pertur- 
bations because these are the ones that couple to the inflaton fluctuations in the linear theory. We 
note that this is not so if we were to consider inflaton fluctuations beyond the linear approxima- 
tion, then tensorial and vectorial metric perturbations would also be driven. The perturbed metric 
9ab = 9ab + ^ab Can be written in the longitudinal gauge as, 

ds'^ = a2(??)[-(l + 2$(x))dr/^ + (1 - 2^{x))Sijdx'dx^], (7.3) 

where the scalar metric perturbations and ^(x) correspond to Bardeen's gauge invariant 

variables [243]. 

7.2 Einstein-Langevin equation for scalar metric perturbations 

The Einstein-Langevin equation as described in section 3 is gauge invariant, and thus we can work 
in a desired gauge and then extract the gauge invariant quantities. The Einstein-Langevin equation 
(3.14) reads now: 

g1? - 87rG(fi°)) + Gi'^ih) - 87TG{f;^l\h)) = SttGU, (7.4) 

where the two first terms cancel, that is G^^^ — 87rG(T^^|j''*) = 0, as the background metric satisfies 
the semiclassical Einstein equations. Here the subscripts (0) and (1) refer to functions in the 
background metric gab and linear in the metric perturbation hab, respectively. The stress tensor 
operator T^b for the minimally coupled inflaton field in the perturbed metric is: 

fab = ya4>^b4> + ^5a(.(Vc0V^(^ + m^^). (7.5) 

Using the decomposition of the scalar field into its homogeneous and inhomogcncous part, sec 
Eq. (7.2), and the metric gab into its homogeneous background gab and its perturbation hab, the 
renormalized expectation value for the stress-energy tensor operator can be written as 

{f^M) = {fab[g])H + {fab[9])<l>^ + {f^b[9])ff, (7.6) 

where the subindices indicate the degree of dependence on the homogeneous field (j) and its pertur- 
bation (p. The first term in this equation depends only on the homogeneous field and it is given by 
the classical expression. The second term is proportional to ((^[y]) which is not zero because the 
field dynamics is considered on the perturbed spacetime, i.e., this term includes the coupling of 
the field with hab and may be obtained from the expectation value of the linearized Klein-Gordon 
equation, 

-m^)^ = 0. (7.7) 

The last term in Eq. (7.6) corresponds to the expectation value to the stress tensor for a free scalar 
field on the spacetime of the perturbed metric. 

After using the previous decomposition, the noise kernel Nahcd[9]x,y) defined in Eq. (3.11) can 
be written as 

{{iab[9;x),icd[9\y)]) = {{iah[9\x),i^d[9]v)]){4>^)2 

x),icd[9\ y)]){^^f , (7.8) 
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where we have used the fact that {(p) = = {0(f"f>) for Gaussian states on the background geometry. 
We consider the vacuum state to be the Euchdean vacuum which is preferred in the de Sitter 
background, and this state is Gaussian. In the above equation the first term is quadratic in (p 
whereas the second one is quartic, both contributions to the noise kernel are separately conserved 
since both (p{r}) and (p satisfy the Klein-Gordon field equations on the background spacetime. 
Consequently, the two terms can be considered separately. On the other hand if one treats as a 
small perturbation the second term in (7.8) is of lower order than the first and may be consistently 
neglected, this corresponds to neglecting the last term of Eq. (7.6). The stress tensor fluctuations 
due to a term of that kind were considered in ref . [244] . 

We can now write down the Einstein-Langevin equations (7.4) to linear order in the inflaton 
fluctuations. It is easy to check [112] that the space-space components coming from the stress tensor 
expectation value terms and the stochastic tensor are diagonal, i.e. {Tij) = = for i / j. This, 
in turn, implies that the two functions characterizing the scalar metric perturbations are equal: 
$ = ^ in agreement with ref. [239] . The equation for $ can be obtained from the Oz-component of 
the Einstein-Langevin equation, which in Fourier space reads 

2iki{n^k + ^'k) = 8TTG{Coi)k, (7.9) 

where ki is the comoving momentum component associated to the comoving coordinate x*, and we 
have used the definition ^k{v) = / d^xe^p{—ik ■ x)^{i],x). Here primes denote derivatives with 
respect to the conformal time rj and TC = a' /a. A nonlocal term of dissipative character which 
comes from the second term in Eq. (7.6) should also appear on the left hand side of Eq. (7.9), 
but we have neglected it to simplify the forthcoming expressions. Its inclusion does not change 
the large scale spectrum in an essential way [112]. Note, however, that the equivalence of the 
stochastic approach to linear order in ip and the usual linear cosmological perturbations approach 
is independent of that approximation [112]. To solve Eq. (7.9), whose left-hand side comes from 
the linearized Einstein tensor for the perturbed metric [239], we need the retarded propagator for 
the gravitational potential 

where / is a homogeneous solution of Eq. (7.9) related to the initial conditions chosen and mp = 
1/G. For instance, if we take f{r],r]') = —6{r]o — r)')a{r]')/a{r]) the solution would correspond to 
"turning on" the stochastic source at rjo. With the solution of the Einstein-Langevin equation (7.9) 
for the scalar metric pcrtiirbations we are in the position to compute the two-point correlation 
functions for these perturbations. 

7.3 Correlation functions for scalar metric perturbations 

The two-point correlation function for the scalar metric perturbations induced by the inflaton 
fluctuations is thus given by 

(<l>fc(??)$fc'(r?'))s = {2nfS{k + k') 
X dr]i dr]2Gk{r],r]i)Gk'{r] ,m){{Coi)k{m){Coi)k'{m))s- (7.11) 
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Here two-point correlation function for the stochastic source, which is connected to the stress-energy 
tensor fluctuations through the noise kernel is given by, 

= lkikici>'{m)(t>'{m)G'i\m,m), (7.i2) 

where 0^^^(771,772) = ({<^fc('7i)! <^-ik(^2)}) is the fe-mode Hadamard function for a free minimally 
coupled scalar field in the appropriate vacuum state on the Priedmann-Robertson- Walker back- 
ground. 

In practice, to make the explicit computation of the Hadamard function we will assume that the 
field state is in the Euclidean vacuum and the background spacetime is de Sitter. Furthermore we 
will compute the Hadamard function for a massless field, and will make a perturbative expansion 
in terms of the dimensionless parameter m/mp. Thus we consider 

G^k\^l^m) = (0|{yfc(77l),y_fc(772)}|0) =27^(7Xfc(77lK(772)), 

with yk{rii) = a{r])(pk{ni) = akUkirj) + aLk''^-k{v) and where 

«fc = (2fe)-i/V'=''(l-V^), 

are the positive frequency fc-mode for a massless minimally coupled scalar field on a de Sitter 
background, which define the Euclidean vacuum state: afc|0) =0 [32]. 

The assumption of a massless field for the computation of the Hadamard function is made 
because massless modes in de Sitter are much simpler to deal with than massive modes. Wc can 
see that this is, however, a reasonable approximation as follows. For a given mode the m = 
approximation is reasonable when its wavelength A is shorter that the Compton wavelength, Ac = 
l/?7i. In our case we have a very small mass m and the horizon size H~^, where H is the Hubble 
constant H = a/a (here a{t) with t the physical time dt = adrj) satisfies that < Ac- Thus, for 
modes inside the horizon A < A^; and rn = is a reasonable approximation. Outside the horizon 
massive modes decay in amplitude as ^ cxp^—rn^t/H) whereas massless modes remain constant, 
thus when modes leave the horizon the approximation will eventually break down. However, we 
only need to ensure that the approximation is still valid after 60 e-folds, i.e. Ht ^ 60, but this is 
the case since 60 < given that m ~ 10~^r7Zp, and m ^ H as m. most inflationary models 
[237, 238]. 

The background geometry is not exactly that of de Sitter spacetime, for which 0(77) = —{Hri)~^ 
with —00 < T] < 0. One can expand in terms of the "slow-roll" parameters and assume that 
to first order <p{t) ~ m'p{m/mp), where t is the physical time. The correlation function for the 
metric perturbation (7.11) can then be easily computed; see Ref. [240, 112] for details. The final 
result, however, is very weakly dependent on the initial conditions as one may understand from 
the fact that the accelerated expansion of de quasi-de Sitter spacetime during inflation erases the 
information about the initial conditions. Thus one may take the initial time to be r/o = —00 and 
obtain to lowest order in m/mp the expression 

{^k{ri)^k'{'n'))s-^T^^ (— )\-H'2T^fS{k + k') COS k{ri-ri'). (7.13) 
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Prom this result two main conclusions are derived. First, the prediction of an almost Harrison- 
Zel'dovich scale-invariant spectrum for large scales, i.e. small vahics of k. Second, since the 
correlation function is of order of {m/mp)^ a severe bound to the mass m is imposed by the grav- 
itational fluctuations derived from the small values of the Cosmic Microwave Background (CMB) 
anisotropics detected by COBE. This bound is of the order of {m/mp) ^ 10~^ [242, 239]. 

We should now comment on some differences with those works in Ref. [108, 109, 110, 111] 
which used a self-interacting scalar field or a scalar field interacting nonlinearly with other fields. 
In those works an important relaxation of the ratio m/nip was found. The long wavelength modes 
of the inflaton field were regarded as an open system in an environment made out of the shorter 
wavelength modes. Then, Langevin type equations were used to compute the correlations of the 
long wavelength modes driven by the fiuctuations of the shorter wavelength modes. In order to 
get a significant relaxation on the above ratio, however, one had to assume that the correlations of 
the free long wavelength modes, which correspond to the dispersion of the system initial state, had 
to be very small. Otherwise they dominate by several orders of magnitude those fluctuations that 
come from the noise of the environment. This would require a great amount of fine-tuning for the 
initial quantum state of each mode [112]. We should remark that in the model discussed here there 
is no environment for the inflaton fluctuations. The inflaton fluctuations, however, are responsible 
for the noise that induce the metric perturbations. 

7.4 Discussion 

One important advantage of the Einstein-Langevin approach to the gravitational fluctuations in 
inflaton over the approach based on the quantization of the linear perturbations of both the met- 
ric and the inflaton field [239], is that an exact treatment of the infiaton quantum fluctuations is 
possible. This leads to corrections to the almost scale invariant spectrum for scalar metric pertur- 
bations at large scales, and has implications for the spectrum of the cosmic microwave background 
anisotropics. However, in the standard inflationary models these corrections are subdominant. 
Furthermore when the full non linear effect of the quantum field is considered, tensorial metric per- 
turbations are also induced by the inflaton fluctuations. An estimation of this effect, presumably 
subdominant over the free tensorial fiuctuations, has not been performed. 

We should remark that although the gravitational fluctuations are here assumed to be classical, 
the correlation functions obtained correspond to the expectation values of the symmetrized quantum 
metric perturbations [187, 112]. This means that even in the absence of decoherence the fluctuations 
predicted by the Einstein-Langevin equation, whose solutions do not describe the actual dynamics 
of the gravitational field any longer, still give the correct symmetrized quantum two-point functions. 

Another important advantage of the stochastic gravity approach is that one may also compute 
the gravitational fluctuations in inflationary models which are not driven by an inflaton field, 
such as Starobinsky inflation which is driven by the trace anomaly due to conformally coupled 
quantum fields. In fact, Einstein semiclassical equation (3.7) for a massless quantum field which 
is conformally coupled to the gravitational field admits an inflationary solution which is almost 
de Sitter initially and ends up in a matter-dominated like regime [113, 114]. In these models the 
standard approach based on the quantization of the gravitational and the matter fields to linear 
order cannot be used. This is because the calculation of the metric perturbations correspond to 
having only the last term in the noise kernel in Eq. (7.8), since there is no homogeneous field (f){ri) 
as the expectation value (0) = 0, and linearization becomes trivial. 
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In the trace anomaly induced inflation framework Hawking et al. [115] were able to compute the 
two-point quantum correlation function for scalar and tcnsorial metric perturbations in a spatially 
closed de Sitter universe, making use of the anti-de Sitter conformal field theory correspondence. 
They find that short scale metric perturbations are strongly suppressed by the conformal matter 
fields. This is similar to what we obtained in section 6 for the induced metric fluctuations in 
Minkowski spacetime. In the stochastic gravity context, the noise kernel in a spatially closed de 
Sitter background was derived in Ref. [244]. But in a spatially fiat arbitrary Friedmann- Robertson- 
Walker model the Einstein-Langevin equation describing the metric perturbations was first obtained 
in Ref. [89], see also [4]. The two-point correlation functions for the metric perturbations can be 
derived from its solutions, but this is work still in progress. 

8 Black Hole Backreaction 

As another illustration of the application of stochastic gravity we consider fiuctuations and backre- 
action in black hole spacetimes. The celebrated Hawking effect of particle creation from black holes 
is constructed from a quantum field theory in curved spacetime framework. The oft-mentioned 
'black hole evaporation' referring to the reduction of the mass of a black hole due to particle cre- 
ation must entail backreaction considerations. Backreaction of Hawking radiation [245, 246, 190, 
191, 192, 247, 248, 249] could alter the evolution of the background spacetime and change the nature 
of its end state, more drastically so for Planck size black holes. Because of the higher symmetry in 
cosmological spacetimes, backreaction studies of processes therein have progressed further than the 
corresponding black hole problems, which to a large degree is still concerned with finding the right 
approximations for the regularized energy momentum tensor [250, 251, 252, 253, 254, 255, 256] 
for even the simplest spacetimes such as the spherically symmetric family including the impor- 
tant Schwarzschild metric (for a summary of the cosmological backreaction problem treated in 
the stochastic gravity theory, see [4]). The latest important work is that of Hiscock, Larson and 
Anderson [256] on backreaction in the interior of a black hole, where one can find a concise sum- 
mary of earlier work. To name a few of the important landmarks in this endeavor (this is adopted 
from [256]), Howard and Candelas [257, 258] have computed the stress-energy of a conformally 
invariant scalar field in the Schwarzschild geometry. Jensen and Ottewill [259] have computed the 
vacuum stress-energy of a massless vector field in Schwarzschild. Approximation methods have 
been developed by Page, Brown, and Ottewill [260, 261, 262] for conformally invariant fields in 
Schwarzschild spacetime, Prolov and Zel'nikov [263] for conformally invariant fields in a general 
static spacetime, Anderson, Hiscock and Samuel [253, 254] for massless arbitrarily coupled scalar 
fields in a general static spherically symmetric spacetime. Furthermore the DeWitt-Schwinger 
approximation has been derived by Frolov and Zel'nikov [264, 265] for massive fields in Kerr 
spacetime, Anderson Hiscock and Samuel [253, 254] for a general (arbitrary curvature coupling 
and mass) scalar field in a general static spherically symmetric spacetime and have applied their 
method to the Reissner-Nordstrom geometry [255]. Though arduous and demanding, the effort 
continues on because of the importance of backreaction effects of Hawking radiation in black holes. 
They are expected to address some of the most basic issues such as black hole thermodynamics 
[147, 148, 149, 266, 267, 268, 269, 270, 271, 272, 273, 274, 275, 276, 277, 278, 279] and the black 
hole end-state and information loss puzzles [280]. 

Here we wish to address the black hole backreaction problem with new insights provided by 
stochastic semiclassical gravity. (For the latest developments see reviews, e.g., [86, 2, 3, 4]). It 
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is not our intention to seek better approximations for the regularized energy momentum tensor, 
but to point out new ingredients lacking in the existing framework based on semiclassical gravity. 
In particular one needs to consider both the dissipation and the fluctuations aspects in the back 
reaction of particle creation or vacuum polarization. 

In a short note [120] Raval, Sinha and one of us discussed the formulation of the problem 
in this new light, commented on some shortcomings of existing works, and sketched the strategy 
[123] behind the stochastic gravity theory approach to the black hole fluctuations and backreaction 
problem. Here we follow their treatment with focus on the class of quasi-static black holes. 

Prom the new perspective provided by statistical field theory and stochastic gravity, it is not dif- 
ficult to postulate that backreaction effect is the manifestation of a fluctuation-dissipation relation 
[281, 282, 283, 284, 285, 286]. This was first conjectured by Candelas and Sciama [116, 117, 118] for 
a dynamic Kerr black hole emitting Hawking radiation, and Mottola [119] for a static black hole (in 
a box) in quasi-equilibrium with its radiation via linear response theory [287, 288, 289, 290, 291]. 
While the fluctuation-dissipation relation in a linear response theory captures the response of the 
system (e.g., dissipation of the black hole) to the environment (in these cases the matter field) 
linear response theory (in the way it is commonly presented in statistical thermodynamics) cannot 
provide a full description of self-consistent backreaction on at least two counts: First, because it 
is usually based on the assumption of a specified background spacctime (static in this case) and 
state (thermal) of the matter field(s)(e.g., [119]). The spacetime and the state of matter should 
be determined in a self-consistent manner by their dynamics and mutual influence. Second, the 
fluctuation part represented by the noise kernel is amiss (e.g., [174, 175]) This is also a problem in 
the fluctuation-dissipation relation proposed by Candelas and Sciama [116, 117, 118] (see below). 
As will be shown in an explicit example later, the back reaction is intrinsically a dynamic process. 
The Einstein-Langevin equation in stochastic gravity overcomes both of these deficiencies. 

For Candelas and Sciama [116, 117, 118], the classical formula they showed relating the dis- 
sipation in area linearly to the squared absolute value of the shear amplitude is suggestive of a 
fluctuation-dissipation relation. When the gravitational perturbations are quantized (they choose 
the quantum state to be the Unruh vacuum) they argue that it approximates a flux of radiation 
from the hole at large radii. Thus the dissipation in area due to the Hawking flux of gravitational 
radiation is allegedly related to the quantum fluctuations of gravitons. The criticism in Rcf. [120] 
is that their's is not an fluctuation-dissipation relation in the truly statistical mechanical sense 
because it does not relate dissipation of a certain quantity (in this case, horizon area) to the fluctu- 
ations of the same quantity. To do so would require one to compute the two point function of the 
area, which, being a four-point function of the graviton field, is related to a two-point function of 
the stress tensor. The stress tensor is the true "generalized force" acting on the spacetime via the 
equations of motion, and the dissipation in the metric must eventually be related to the fluctuations 
of this generalized force for the relation to qualify as an fluctuation-dissipation relation. 

From this reasoning, we see that the stress energy bi-tensor and its vacuum expectation value 
known as the noise kernel, arc the new ingredients in backreaction considerations. But these are 
exactly the centerpiece in stochastic gravity. Therefore the correct framework to address semi- 
classical backreaction problems is stochastic gravity theory, where fluctuations and dissipation 
arc the equally essential components. The noise kernel for quantum fields in Minkowski and de 
Sitter spacetime has been carried out by Martin, Roura and Verdaguer [5, 94, 112], for thermal 
fields in black hole spacetime and scalar fields in general spacetimes by Campos, Hu and Phillips 
[121, 122, 7, 8, 167]. Earlier, for cosmological backreaction problems Hu and Sinha [85] derived 
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a generalized expression relating dissipation (of anisotropy in Bianchi Type I universes) and fluc- 
tuations (measured by particle numbers created in neighboring histories). This example shows 
that one can understand the backreaction of particle creation as a manifestation of a (generalized) 
fluctuation-dissipation relation. 

As an illustration of the application of stochastic gravity theory we outline the steps in a 
black hole backreaction calculation, focusing on the manageable quasi-static class. We adopt the 
Hartle-Hawking picture [292] where the black hole is bathed eternally - actually in quasi-thermal 
equilibrium - in the Hawking radiance it emits. It is described here by a massless scalar quantum 
field at the Hawking temperature. As is well-known, this quasi-equilibrium condition is possible 
only if the black hole is enclosed in a box of size suitably larger than the event horizon. We 
can divide our consideration into the far field case and the near horizon case. Campos and Hu 
[121, 122] have treated a relativistic thermal plasma in a weak gravitational field. Since the far 
field limit of a Schwarzschild metric is just the perturbed Minkowski spacetimc, one can perform a 
perturbation expansion off hot flat space using the thermal Green functions [293]. Strictly speaking 
the location of the box holding the black hole in equilibrium with its thermal radiation is as far 
as one can go, thus the metric may not reach the perturbed Minkowski form. But one can also 
put the black hole and its radiation in an anti-de Sitter space [294], which contains such a region. 
Hot flat space has been studied before for various purposes. See e.g., [295, 296, 297, 298, 299]. 
Campos and Hu derived a stochastic CTP effective action and from it an equation of motion, the 
Einstein Langevin equation, for the dynamical effect of a scalar quantum field on a background 
spacetime. To perform calculations leading to the Einstein-Langevin equation one needs to begin 
with a self-consistent solution of the semiclassical Einstein equation for the thermal field and the 
perturbed background spacetime. For a black hole background, a semiclassical gravity solution is 
provided by York [190, 191, 192]. For a Robertson- Walker background with thermal fields it is given 
by Hu [300]. Recently Sinha, Raval and Hu [123] outlined a strategy for treating the near horizon 
case, following the same scheme of Campos and Hu. In both cases two new terms appear which 
are absent in semiclassical gravity considerations: a nonlocal dissipation and a (generally colored) 
noise kernel. When one takes the noise average one recovers York's [190, 191, 192] semiclassical 
equations for radially perturbed quasi-static black holes. For the near horizon case one cannot 
obtain the full details yet, because the Green function for a scalar field in the Schwarzschild metric 
comes only in an approximate form (e.g. Page approximation [260]), which, though reasonably 
accurate for the stress tensor, fails poorly for the noise kernel [8, 167]. In addition a formula is 
derived in [123] expressing the CTP effective action in terms of the Bogolyubov coefficients. Since 
it measures not only the number of particles created, but also the difference of particle creation in 
alternative histories, this provides a useful avenue to explore the wider set of issues in black hole 
physics related to noise and fluctuations. 

Since backreaction calculations in semiclassical gravity has been under study for a much longer 
time than in stochastic gravity we will concentrate on explaining how the new stochastic features 
arise from the framework of semiclassical gravity, i.e., noise and fluctuations and their conse- 
quences. Technically the goal is to obtain an influence action for this model of a black hole coupled 
to a scalar field and to derive an Einstein-Langevin equation from it. As a by-product, from the 
fluctuation-dissipation relation, one can derive the vacuum susceptibility function and the isother- 
mal compressibility function for black holes, two quantities of fundamental interest in characterizing 
the nonequilibrium thermodynamic properties of black holes. 
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8.1 The Model 



In this model the black hole spacetime is described by a spherically symmetric static metric with line 
element of the following general form written in advanced time Eddington-Finkelstein coordinates 

ds^ = g^ydx^dx" = -e^^ - dv"^ + 2e^^dvdr + d^'^, (8.1) 

where if) = il){r) and m = m{r) , v = t + r + 2Mln — 1) and dO^ is the line element on the two 
sphere. Hawking radiation is described by a massless, conformally coupled quantum scalar field (f) 
with the classical action 

Sm[<f^, g^u] = -\j d^x^[gi'''d^<t>dA + ^{ri)Rct>\ (8.2) 

where ^(n) = {n is the dimension of spacetime) and R is the curvature scalar of the spacetime 

it lives in. 

Let us consider linear perturbations h^i, off a background Schwarzschild metric g^J 

5M- = 5i°J+V. (8-3) 

with standard line element 

We look for this class of perturbed metrics in the form given by (8.1), (thus restricting our consid- 
eration only to spherically symmetric perturbations): 

~ 1 + ep{r), (8.5) 

and 

m M[l + en{r)], (8.6) 

where = \aTf^; a = |^;A = 90(8^)7r^. Th is the Hawking temperature. This particular 
parametrization of the perturbation is chosen following York's [190, 191, 192] notation. Thus the 
only non-zero components of h^i, are 

2M, , , 2Mea(r)\ , , 

Kv = -[i^ 2ep r + ^ , 8.7 

and 

Kr = ep{r). (8.8) 

So this represents a metric with small static and radial perturbations about a Schwarzschild black 
hole. The initial quantum state of the scalar field is taken to be the Hartle Hawking vacuum, which 
is essentially a thermal state at the Hawking temperature and it represents a black hole in (unstable) 
thermal equilibrium with its own Hawking radiation. In the far field limit, the gravitational field 
is described by a linear perturbation from Minkowski spacetime. In equilibrium the thermal bath 
can be characterized by a relativistic fluid with a four-velocity (time-like normalized vector field) 
M^, and temperature in its own rest frame (3^^. 
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To facilitate later comparisons with our program we briefly recall York's work [190, 191, 192]. 
See also work by Hocliberg and Kephart [247] for a massless vector field, Hochberg, Kephart and 
York [248] for a massless spinor field, and Anderson, Hiscock, Whitesell, and York [249] for a 
quantized massless scalar field with arbitrary coupling to spacetime curvature. York considered the 
semiclassical Einstein equation 



with Gfj,i, ~ G^iJ +6Gfj,i, where G^iJ is the Einstein tensor for the background spacetime. The zeroth 
order solution gives a background metric in empty space, i.e, the Schwarzschild metric. 5Gfj_i, is 
the linear correction to the Einstein tensor in the perturbed metric. The semiclassical Einstein 
equation in this approximation therefore reduces to 



York solved this equation to first order by using the expectation value of the energy momentum 
tensor for a conformally coupled scalar field in the Hartle-Hawking vacuum in the unperturbed 
(Schwarzschild) spacetime on the right hand side and using (8.7) and (8.8) to calculate 5G^i^ on 
the left hand side. Unfortunately, no exact analytical expression is available for the (T^i^) in a 
Schwarzschild metric with the quantum field in the Hartle-Hawking vacuum that goes on the right 
hand side. York therefore uses the approximate expression given by Page [260] which is known 
to give excellent agreement with numerical results. Page's approximate expression for {T^^) was 
constructed using a thermal Feynman Green's function obtained by a conformal transformation of 
a WKB approximated Green's function for an optical Schwarzschild metric. York then solves the 
semiclassical Einstein equation (8.10) self consistently to obtain the corrections to the background 
metric induced by the backreaction encoded in the functions /u(r) and p{r). There was no mention of 
fluctuations or its effects. As we shall see, in the language of the previous section, the semiclassical 
gravity procedure which York followed working at the equation of motion level is equivalent to 
looking at the noise-averaged backreaction effects. 

8.2 CTP Effective Action for the Black Hole 

We first derive the CTP effective action for the model described in the previous section. Using the 
metric (8.4) (and neglecting the surface terms that appear in an integration by parts) we have the 
action for the scalar field written perturbatively as 



(8.9) 



6G^,{g^''\h) = K{T^,). 



(8.10) 



(8.11) 



where the first and second order perturbative operators F^^^ and V^'^^ are given by 





a,-^(n)(i2(2) + i/ii2W)|. 



(8.12) 
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In the above expressions, R^'^^ is the /c-order term in the perturbation hi^i,{x) of the scalar curvature 
R and h/^i, and denote a hnear and a quadratic combination of the perturbation, respectively, 

h = h - l/io(0) 
"■//I/ — "'fiu 2 ^M'^ ' 

V = h^hau-lhh^. + lh^9i^J -Ihaph'^^gj^l (8.13) 

From quantum field theory in curved spacetime considerations discussed above we take the following 
action for the gravitational field: 



3R^,^fs{x)R^''^f'{x)- 



1-360 



R^{x)}. (8.14) 



The first term is the classical Einstein-Hilbert action and the second term is the counterterm in 
four dimensions used to renormalize the divergent effective action. In this action = 167rGAr, 
a = (28807r2)-i and /X is an arbitrary mass scale. 

We are interested in computing the CTP effective action (8.11) for the matter action and when 
the field (p is initially in the Hartle-Hawking vacuum. This is equivalent to saying that the initial 
state of the field is described by a thermal density matrix at a finite temperature T = Th ■ The CTP 
effective action at finite temperature T = \/f3 for this model is given by (for details see [121, 122]) 

sLeffHu] = S,AhU - S,[h;J-'-Tr{\nG^Jh%]}, (8.15) 



where lb denote the forward and backward time path of the CTP formalism and G^^i,[h^A is the 

complete 2x2 matrix propagator (a and 6 take it values: and G correspond to 

the Feynman, Wightman and Schwinger Greens functions respectively) with thermal boundary 

conditions for the differential operator + • • •). The actual form of G^^ 

cannot be explicitly given. However, it is easy to obtain a perturbative expansion in terms of 

ik') . ik') (k) 

V^f^ , the A;-order matrix version of the complete differential operator defined by = ±V^ 

and = 0, and G^^, the thermal matrix propagator for a massless scalar field in Schwarzschild 
spacetime . To second order reads, 

Gf, = Gf, - Gf^K^Gi - GlV^^Gi + GlV^^Giv^fG% + • • • (8.16) 

Expanding the logarithm and dropping one term independent of the perturbation /i^^ {x) , the CTP 
effective action may be perturbatively written as, 

^\Tr[vi'^Gi^-V^'^Gt_+vi^^GU-V^yGU 

-jrr[y«G^+vi^)G^++yi^)G^_y«G^_ 

-2v!^^Gl_V^^Gt^\ (8.17) 

In computing the traces, some terms containing divergences are canceled using counterterms intro- 
duced in the classical gravitational action after dimensional regularization. 
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8.3 Near Flat Case 

At this point we divide our considerations into two cases. In the far field limit h^i^ represent 

perturbations about flat space, i.e., g^J = rjfii,. The exact "unperturbed" thermal propagators 

for scalar fields are known, i.e., the Euclidean propagator with periodicity (3. Using the Fourier 
transformed form (those quantities are denoted with a tilde) of the thermal propagators G^f^{k) , 

the trace terms of the form Tr[Va^^ Gf^^V^^^ Gl^g] can be written as [121, 122], 
TrM'^GiX'^G^J = I «V hl,{x)hip{x') 

-^^''''^^^^^^'^nik + <l)&rsi^^^^^^ (8-18) 

where the tensor T^'^'"/^ (q, k) is defined in [121, 122] after an expansion in terms of a basis of 
14 tensors [296, 297]. In particular, the last trace of (8.17) may be split in two different kernels 
W''^'^{x - x') and B^"'^'^f^{x - x'), 

'-Tr[vi'^Gi_vPG^^] = 

- J d^xd^x' /i+ (x)/t;^(x')[D^'^'°^(x - x') + m^^^'^^{x - x')]. (8.19) 
One can express the Fourier transforms of these kernels, respectively, as 

N'^''"''(A;)=7r2y ^ {e{^k" + qO)e{-qO)+e{-k'>-q'>)e{q") 

+npm)+np{\k° + q°\) 

+2npm)np{\k° + q"\)] 5{q^)5[{k + qf]T^^'''''^{q, k), (8.20) 

B"''''\k)=-in' I ^ {e{k^ + qO)e{-q'^)-e{-k'^-q")e{q°) 
+S9{k" + q")np{\q''\) 

-sg{q'')np{r + g1)} 5{q^mk + g)2]T'^-'«/^(g, k). (8.21) 

Using the property Ti^^'<^P{q, k) = T'^'"''^'^ {-q, -k), it is easy to see that the kernel N^''''°'^{x - x') is 
symmetric and D^'^'"'^(.t — x') antisymmetric in their arguments; that is, N'*'''"''(a;) = 'N^^''^^{—x) 
and D^'^'°^(a;) = -D^'^'"^(-x). 

The physical meanings of these kernels can be extracted if we write the renormalized OTP 
effective action at finite temperature (8.17) in an influence functional form [301, 302, 303, 304]. N, 
the imaginary part of the CTP effective action can be identified with the noise kernel and D, the 
antisymmetric piece of the real part, with the dissipation kernel. Campos and Hu [121, 122] have 
shown that these kernels identified as such indeed satisfy a thermal fluctuation-dissipation relation. 

If we denote the difference and the sum of the perturbations /i^j, deflned along each branch C± 
of the complex time path of integration C by [/i^jy] = /i^,^ — and {h^^} = h^^ + h'^^, respectively, 
the influence functional form of the thermal CTP effective action may be written to second order 
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in ha,y as, 



^-ff^^J'^l ^ 2(16^ / [V](^)Lf:r^(^ - x'){hap}{x') 



+^ j d^x d'^x' [V](^)H''"'"^(x - x'){Kf)}{x') 
-\ J d^x d\' [h^,]{xW'^^{x - x'){Kp}{x') 
+ j d'^x d\' [V](^)N'^"'"^(x - x')[Kp\{x'). 



The first fine is the Einstein-Hilbert action to second order in the perturbation h^^{x). L^^^' 



(o) 



(8.22) 



is a symmetric kernel (i.e. h'f^-!'^^(x) = L/'^'"'^(— x)). In the near flat case its Fourier transform is 



given by 



Ho) 



(k) 



1 



-ik'Tf '"^(^, k) + 2A;'Tr"^(g, A;) 
+Tp°''(g, k)-2Tl^'''^{q,k) 



(8.23) 



The fourteen elements of the tensor basis T^^'° (g, k) (i = 1, • • • , 14) are defined in [296, 297]. The 
second is a local term linear in h^^{x). Only when far away from the hole that it takes the form of 
the stress tensor of massless scalar particles at temperature which has the form of a perfect 
fluid stress-energy tensor 



TT 



rpfJ-f 

" 30/34 



?xV + ^(r?'^^ + uV) 



(8.24) 



where u'^ is the four-velocity of the plasma and the factor is the familiar thermal energy 
density for massless scalar particles at temperature In the far field limit, taking into account 
the four-velocity of the fluid, a manifestly Lorentz-covariant approach to thermal fleld theory 
may be used [305]. However, in order to simplify the involved tensorial structure we work in the 
co-moving coordinate system of the fluid where = (1,0,0,0). In the third line, the Fourier 
transform of the symmetric kernel 'R^'^''^^(x) can be expressed as 



ii'"'^''^{k) 



a 



k^ f 1 



^ In ^Q'^'^>°/3(A;) + ^Q'^'''"'^(A;) 



+^A{-'^'i'''''^^^k)-2T:^^''^^u,k) 

+Tr"^(^/,yt) + 2Tf'"^(n,fe)} 

+^ {fc'Tr"''(^, k) - 2k'Tr-^{u, k) 

-T^'^'"^(n, k) + 2T:'il^''\u, k)} 
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+5[{k+q)^]np{\k"+q"\)V X T'^-'°/3(g, fc), (8.25) 



where is a simple redefinition of the rcnornialization parameter fi given by // = iJ,ex.p[^ + 2 ln47r — 
57), and the tensors Q^'^'"^{k) and Q,'^'^'"'^ {k) are defined, respectively, by 

QM.,a/3(^) = ^ {Tr"^(9, k) - ^Tr"^(g, fc) + ^.Tl^'^'^iq, A:)} 

_[l_360(e-^)2]{Tr"''(g,fe) + ^T^2''"^(g,fe)-pT^3''"''(g,fc)}^ (8.26) 

Q'''''"^(fc)=[l + 576(e - If - 60(^ - 1)(1 - 36^')] {Tr"^(g, fc) 

+^T^2'"^('?, fe) - ^T^'ir^Q, fc)} • (8.27) 

In the above and subsequent equations, we denote the coupling parameter in four dimensions ^(4) by 
^ and consequently ^' means d^{n)/dn evaluated at n = 4. h'*^'" (fc) is the complete contribution of 
a free massless quantum scalar field to the thermal graviton polarization tensor[296, 297, 298, 299] 
and it is responsible for the instabilities found in flat spacetime at finite temperature [295, 296, 
297, 298, 299]. Note that the addition of the contribution of other kinds of matter fields to the 
effective action, even graviton contributions, does not change the tensor structure of these kernels 
and only the overall factors are different to leading order [296, 297]. Eq. (8.25) rcfiects the fact 

that the kernel H'"''"'^(A;) has thermal as well as non-thermal contributions. Note that it reduces 
to the first term in the zero temperature limit (/3 — ^ 00) 

H'^-."/3(fc) ^ _^|Ln^Q''-'"/3(fe) + ^Q'^^'"^(A;)|. (8.28) 

and at high temperatures the leading term {(3~*) may be written as 

H"'^'"''(fc) - ^,J>^{r)^^^n,K), (8.29) 

where we have introduced the dimensionless external momentum K^^ = k^^ /\k\ = (r, k). The Hj(r) 
coefficients were first given in [296, 297] and generalized to the next-to-leading order (/?~^) in 
[298, 299]. (They are given with the MTW sign convention in [121, 122].) 

Finally, as defined above, N^'''"'^(x) is the noise kernel representing the random fiuctuations of 
the thermal radiance and D'^'^'°'^{x) is the dissipation kernel, describing the dissipation of energy 
of the gravitational field. 

8.4 Near Horizon Case 

In this case, since the perturbation is taken around the Schwarzschild spacetime, exact expressions 
for the corresponding unperturbed propagators Gf^[/iui/] ^ot known. Therefore apart from the 
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approximation of computing the CTP effective action to certain order in perturbation theory, an 
appropriate approximation scheme for the unperturbed Green's functions is also required. This 
feature manifested itself in York's calculation of backreaction as well, where, in writing the {T^j^y) 
on the right hand side of the semiclassical Einstein equation in the unperturbed Schwarzschild 
metric, he had to use an approximate expression for [T^y) in the Schwarzschild metric given by 
Page [260]. The additional complication here is that while to obtain {T^y) as in York's calculation, 
the knowledge of only the thermal Feynman Green's function is required, to calculate the CTP 
effective action one needs the knowledge of the full matrix propagator, which involves the Feynman, 
Schwinger and Wightman functions. 

It is indeed possible to construct the full thermal matrix propagator G^^^h^y] based on Page's 
approximate Feynman Green's function by using identities relating the Feynman Green's function 
with the other Green's functions with different boundary conditions. One can then proceed to 
explicitly compute a CTP effective action and hence the influence functional based on this approx- 
imation. However, we desist from delving into such a calculation for the following reason. Our 
main interest in performing such a calculation is to identify and analyze the noise term which is 
the new ingredient in the backreaction. We have mentioned that the noise term gives a stochastic 
contribution to the Einstein-Langevin equation (3.14). We had also stated that this term is 
related to the variance of fluctuations in T^y, i.e, schematically, to (T^y). However, a calculation 
of {T^y) in the Hartlc-Hawking state in a Schwarzschild background using the Page approximation 
was performed by Phillips and Hu [7, 8, 167] and it was shown that though the approximation is 
excellent as far as {T^y) is concerned, it gives unacceptably large errors for {T^y) at the horizon. 
In fact, similar errors will be propagated in the non-local dissipation term as well, because both 
terms originate from the same source, that is, they come from the last trace term in (8.17) which 
contains terms quadratic in the Green's function. However, the Influence Functional or CTP for- 
malism itself does not depend on the nature of the approximation, so we will attempt to exhibit 
the general structure of the calculation without resorting to a specific form for the Greens function 
and conjecture on what is to be expected. A more accurate computation can be performed using 
this formal structure once a better approximation becomes available. 

The general structure of the CTP effective action arising from the calculation of the traces in 
equation (8.17) remains the same. But to write down explicit expressions for the non-local kernels 
one requires the input of the explicit form of G^^^lh^y] in the Schwarzschild metric, which is not 
available in closed form. We can make some general observations about the terms in there. The 
first line containing L does not have an explicit Fourier representation as given in the far field 
case, neither will T^^-^ in the second line representing the zeroth order contribution to have 
a perfect fluid form. The third and fourth terms containing the remaining quadratic component 
of the real part of the effective action will not have any simple or even complicated analytic form. 
The symmetry properties of the kernels H^^'^'°'l^{x,x') and D'^'^'°'l^{x,x') remain intact, i.e., they are 
respectively even and odd in x, x'. The last term in the CTP effective action gives the imaginary 
part of the effective action and the kernel N{x, x') is symmetric. 

Continuing our general observations from this CTP effective action, using the connection be- 
tween this thermal CTP effective action to the influence functional [18, 87] via an equation in the 
schematic form (4.3). We see that the nonlocal imaginary term containing the kernel N^"'°'^(x,x') 
is responsible for the generation of the stochastic noise term in the Einstein-Langevin equation 
and the real non-local term containing kernel D'^^''^f^{x,x') is responsible for the non-local dissi- 
pation term. To derive the Einstein-Langevin equation we first construct the stochastic effective 
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action (4.13). We then derive the equation of motion, as shown earHer in (4.15), by taking its 
functional derivative with respect to [hfj_y\ and equating it to zero. With the identification of noise 
and dissipation kernels, one can write down a linear, non-local relation of the form, 

N{t -t')= j d{s - s')K{t -t',s- s')7(s - s'), (8.30) 

where D{t,t') = —dtrj(t,t'). This is the general functional form of a fluctuation-dissipation re- 
lation and K{t,s) is called the fluctuation-dissipation kernel [301, 302, 303, 304]. In the present 
context this relation depicts the backreaction of thermal Hawking radiance for a black hole in 
quasi-equilibr ium . 



8.5 Einstein-Langevin equation 

In this section we show how a semiclassical Einstein-Langevin equation can be derived from the 
previous thermal CTP effective action. This equation depicts the stochastic evolution of the per- 
turbations of the black hole under the influence of the fluctuations of the thermal scalar field. 

The influence functional J-ip = exp(i5iF) previously introduced in Eq. (4.2) can be written 
in terms of the the CTP effective action •S'^fj [/i^j^] derived in equation (8.22) using Eq. (4.3). The 
Einstein-Langevin equation follows from taking the functional derivative of the stochastic effective 
action (4.13) with respect to [hfj,^]{x) and imposing [hn^]{x) = 0. This leads to 

/ d'x' Ll^-^x - x')Kp{x') + \ rg; + / d'x' {w^^-f^{x - x') 

_ ^')^ h^p{x') + ^^'^(x) = 0. (8.31) 

where 

{e''{x)e^{x'))j = W^^^^^ix-x'), (8.32) 

In the far field limit this equation should reduce to that obtained by Campos and Hu [121, 122]: 
For gravitational perturbations h^^'^ defined in (8.13) under the harmonic gauge h^^^y = 0, their 
Einstein-Langevin equation is given by 

_D/^^,a/3(^ - x')) hi^{x') + 2i>'''{x)] = 0, (8.33) 

where the tensor Ppa,ai3 is given by 

Ppa,aP = ^ iVpaVaP + Vp/Sritra - VpaVafi) ■ (8-34) 

The expression for Ppa,ai3 in the near horizon limit of course cannot be expressed in such a simple 
form. Note that this differential stochastic equation includes a non-local term responsible for the 
dissipation of the gravitational field and a noise source term which accounts for the fluctuations 
of the quantum field . Note also that this equation in combination with the correlation for the 
stochastic variable (8.32) determine the two-point correlation for the stochastic metric fiuctuations 
{hjj,uix)hap{x')) ^ self-consistently. 
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As we have seen before and here, the Einstein-Langevin equation is a dynamical equation 
governing the dissipative evolution of the gravitational field under the influence of the fluctuations 
of the quantum field, which, in the case of black holes, takes the form of thermal radiance. From 
its form we can see that even for the quasi-static case under study the back reaction of Hawking 
radiation on the black hole spacetime has an innate dynamical nature. 

For the far field case making use of the explicit forms available for the noise and dissipation 
kernels Campos and Hu [121, 122] formally proved the existence of a fluctuation-dissipation relation 
at all temperatures between the quantum fluctuations of the thermal radiance and the dissipation of 
the gravitational field. They also showed the formal equivalence of this method with linear response 
theory for lowest order perturbations of a near-equilibrium system, and how the response functions 
such as the contribution of the quantum scalar field to the thermal graviton polarization tensor 
can be derived. An important quantity not usually obtained in linear response theory, but of equal 
importance, manifest in the CTP stochastic approach is the noise term arising from the quantum 
and statistical fluctuations in the thermal field. The example given in this section shows that the 
back reaction is intrinsically a dynamic process described (at this level of sophistication) by the 
Einstein-Langevin equation. By comparison, traditional linear response theory calculations cannot 
capture the dynamics as fully and thus cannot provide a complete description of the backreaction 
problem. 

8.6 Discussions 

We make a few remarks here and draw some connection with related work on black hole fluctuations. 
8.6.1 Black hole Backreaction 

As remarked earlier, except for the near-fiat case, an analytic form of the Green function is not 
available. Even the Page approximation [260] which gives unexpectedly good results for the stress 
energy tensor has been shown to fail in the fluctuations of the energy density [8, 167]. Thus using 
such an approximation for the noise kernel will give unreliable results for the Einstein-Langevin 
equation. If we confine ourselves to Page's approximation and derive the equation of motion without 
the stochastic term, we expect to recover York's semiclassical Einstein's equation if one retains only 
the zeroth order contribution, i.c, the first two terms in the expression for the CTP effective action 
in Eq. (8.22). Thus, this offers a new route to arrive at York's semiclassical Einstein's equations. 
Not only is it a derivation of York's result from a different point of view, but it also shows how 
his result arises as an appropriate limit of a more complete framework, i.e, it arises when one 
averages over the noise. Another point worth noting is that our treatment will also yield a non- 
local dissipation term arising from the fourth term in equation (8.22) in the CTP effective action 
which is absent in York's treatment. This difference arises primarily due to the difference in the 
way backreaction is treated, at the level of iterative approximations on the equation of motion as 
in York, versus the treatment at the effective action level as pursued here. In York's treatment, the 
Einstein tensor is computed to first order in perturbation theory , while (T^iy) on the right hand 
side of the semiclassical Einstein equation is replaced by the zeroth order term. In the effective 
action treatment the full effective action is computed to second order in perturbation, and hence 
includes the higher order non-local terms. 

The other important conceptual point that comes to light from this approach is that related 
to the Fluctuation-Dissipation Relation. In the quantum Brownian motion analog (e.g., [301, 
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302, 303, 304] and references therein), the dissipation of the energy of the Brownian particle as 
it approaches equihbriTim and the flTictuations at equihbrium are connected by the Fluctuation- 
Dissipation relation. Here the backreaction of quantum fields on black holes also consists of two 
forms - dissipation and fluctuation or noise - corresponding to the real and imaginary parts of 
the influence functional as embodied in the dissipation and noise kernels. A fluctuation-dissipation 
relation relation has been shown to exist for the near flat case by Campos and Hu [121, 122] and we 
anticipate that it should also exist between the noise and dissipation kernels for the general case, as it 
is a categorical relation [301, 302, 303, 304, 86]. Martin and Verdaguer have also proved the existence 
of a fluctuation-dissipation relation when the semiclassical background is a stationary spacetime 
and the quantum field is in thermal equilibrium. Their result was then extended to a conformal 
field in a conformally stationary background [5]. The existence of a fluctuation-dissipation relation 
for the black hole case has been discussed by some authors previously [116, 117, 118, 119]. In [120], 
Hu, Raval and Sinha have described how this approach and results differ from those of previous 
authors. The fluctuation-dissipation relation reveals an interesting connection between black holes 
interacting with quantum fields and non-equilibrium statistical mechanics. Even in its restricted 
quasi-static form, this relation will allow us to study nonequilibrium thermodynamic properties 
of the black hole under the influence of stochastic fluctuations of the energy momentum tensor 
dictated by the noise terms. 

There are limitations of a technical nature in the specific example invoked here. For one we have 
to confine ourselves to small perturbations about a background metric. For another, as mentioned 
above, there is no reliable approximation to the Schwarzschild thermal Green's function to explicitly 
compute the noise and dissipation kernels. This limits our ability to present explicit analytical 
expressions for these kernels. One can try to improve on Page's approximation by retaining terms 
to higher order. A less ambitious first step could be to confine attention to the horizon and using 
approximations that are restricted to near the horizon and work out the Influence Functional in 
this regime. 

Yet another technical limitation of the specific example is the following. Though we have 
allowed for backreaction effects to modify the initial state in the sense that the temperature of the 
Hartle-Hawking state gets affected by the backreaction, we have essentially confined our analysis 
to a Hartle-Hawking thermal state of the field. This analysis does not directly extend to a more 
general class of states, for example to the case where the initial state of the field is in the Unruh 
vacuum. Thus we will not be able to comment on issues of the stability of an isolated radiating 
black hole under the influence of stochastic fluctuations. 

8.6.2 Metric Fluctuations in Black Holes 

In addition to the work described above by Campos, Hu, Raval and Sinha [121, 122, 120, 123] and 
earlier work quoted therein, we mention also some recent work on black hole metric fluctuations 
and their effect on Hawking radiation. For example, Casher et al [306] and Sorkin [100, 101] have 
concentrated on the issue of fluctuations of the horizon induced by a fluctuating metric. Casher 
et al [306] considers the fluctuations of the horizon induced by the "atmosphere" of high angular 
momentum particles near the horizon, while Sorkin [100, 101] calculates fluctuations of the shape 
of the horizon induced by the quantum field fluctuations under a Newtonian approximations. Both 
group of authors come to the conclusion that horizon fluctuations become large at scales much 
larger than the Planck scale (note Ford and Svaiter [98] later presented results contrary to this 
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claim). However, though these works do deal with backreaction, the fluctuations considered do not 
arise as an explicit stochastic noise term as in our treatment. It may be worthwhile exploring the 
horizon fluctuations induced by the stochastic metric in our model and comparing the conclusions 
with the above authors. Barrabes et al [102, 103] have considered the propagation of null rays and 
massless fields in a black hole fluctuating geometry and have shown that the stochastic nature of the 
metric leads to a modified dispersion relation and helps to confront the trans-Planckian frequency 
problem. However, in this case the stochastic noise is put in by hand and does not naturally arise 
from coarse graining as in the quantum open systems approach. It also does not take backreaction 
into account. It will be interesting to explore how a stochastic black hole metric, arising as a 
solution to the Einstein-Langevin equation, hence fully incorporating backreaction, would affect 
the trans-Planckian problem. 

Ford and his collaborators [98, 173, 99] have also explored the issue of metric fluctuations in 
detail and in particular have studied the fluctuations of the black hole horizon induced by metric 
fluctuations. However, the fluctuations they considered are in the context of a fixed background 
and do not relate to the backreaction. 

Another work originating from the same vein of stochastic gravity but not complying with the 
backreaction spirit is that of Hu and Shiokawa [97] , who study effects associated with electromag- 
netic wave propagation in a Robertson- Walker universe and the Schwarzschild spacetime with a 
small amount of given metric stochasticity. They find that time-independent randomness can de- 
crease the total luminosity of Hawking radiation due to multiple scattering of waves outside the 
black hole and gives rise to event horizon fluctuations and fluctuations in the Hawking tempera- 
ture. The stochasticity in a background metric in their work is assumed rather than derived (from 
quantum field fluctuations, as in this work) and so is not in the same spirit of backreaction. But it 
is interesting to compare their results with that of backreaction, so one can begin to get a sense of 
the different sources of stochasticity and their weights (see, e.g., [2] for a list of possible sources of 
stochasticity.) 

In a subsequent paper Shiokawa [130] showed that the scalar and spinor waves in a stochastic 
spacetime behave similarly to the electrons in a disordered system. Viewing this as a quantum 
transport problem, he expressed the conductance and its fluctuations in terms of a nonlinear sigma 
model in the closed time path formalism and showed that the conductance fluctuations are universal, 
independent of the volume of the stochastic region and the amount of stochasticity. This result 
can have signiflcant importance in characterizing the mesoscopic behavior of spacetimes resting 
between the semiclassical and the quantum regimes. 

9 Concluding Remarks 

In the first part of this review on the fundamentals of theory we have given two routes to the 
establishment of stochastic gravity and derived a general (finite) expression for the noise kernel. In 
the second part we gave three applications, the correlation functions of gravitons in a perturbed 
Minkowski metric, structure formation in stochastic gravity theory and the outline of a program 
for the study of black hole fluctuations and backreaction. A central issue which stochastic gravity 
can perhaps best address is the validity of semiclassical gravity as measured by the fluctuations of 
stress energy compared to the mean. We will include a review of this topic in a future update. 

There are a number of ongoing research related to the topics discussed in this review. On the 
theory side, Roura and Verdaguer [169] has recently showed how stochastic gravity can be related 
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to the large limit of quantum metric fluctuations. Given free matter fields interacting with the 
gravitational field, Hartle and Horowitz [181] and Tomboulis [182] have shown that semiclassical 
gravity can be obtained as the leading order large N limit (while keeping N times the gravitational 
coupling constant fixed). It is of interest to find out where in this setting can one place the 
fluctuations of the quantum fields and the metric fluctuations they induce; specifically, whether 
the inclusion of these sources will lead to an Einstein-Langevin equation [87, 88, 85, 89, 92], as it 
was derived historically in other ways, as described in the first part of this review. This is useful 
because it may provide another pathway or angle in connecting semiclassical to quantum gravity 
(a related idea is the kinetic theory approach to quantum gravity described in [133]). 

Theoretically, stochastic gravity is at the frontline of the 'bottom-up' approach to quantum 
gravity [131, 2, 133]. Structurally, as can be seen from the issues discussed and the applications 
given, stochastic gravity has a very rich constituency because it is based on quantum field theory 
and nonequilibrium statistical mechanics in a curved spacetime context. The open systems concepts 
and the closed-timc-path / influence functional methods constitute an extended framework suitable 
for treating the backreaction and fluctuations problems of dynamical spacetimes interacting with 
quantum fields. We have seen it applied to cosmological backreaction problems. It can also be 
applied to treat the backreaction of Hawking radiation in a fully dynamical black hole collapse 
situation, and in so doing enable one to address related issues such as the black hole end state and 
information loss puzzles (e.g., [280, 307] and references therein). The main reason why this program 
has not progressed as swiftly as desired is due more to technical rather than programatic difficulties 
(such as finding reasonable analytic approximations for the Green function or numerical evaluation 
of mode-sums near the black hole horizon). Finally, the multiplex structure of this theory could be 
used to explore new lines of inquiry and launch new programs of research, such as nonequilibrium 
black hole thermodynamics and statistical mechanics. 
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